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ABSTRACT 
This  is  a  report  on  atmospheric  boundary  layer  measurements  made 
by  NPS  personnel  aboard  the  USNS  Kane  in  February  and  March  of  1978. 
The  data  has  been  analyzed  to  obtain  Monin-Obukhov  similarity  para- 
meters and  EM  propagation  properties  using  bulk,  profile  and  turbu- 
lence data.  The  evaporation  duct  height,  Z*,  was  found  to  vary  from 
1.0  meters  to  25  meters  with  a  typical  value  of  16  meters.  A  compari- 
son with  profile  and  turbulence  methods  showed  that  the  bulk  method 
was  accurate  to  within  a  factor  of  two,  although  some  of  the  disagree- 
ment has  been  attributed  to  profile  measurement  inaccuracies.  A  sen- 
sitivity analysis  showed  that  the  relative  humidity  and  air-sea  tem- 
perature difference  are  the  most  critical  measurements  for  bulk  calcu- 
lations of  Z*. 
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I.       Summary 

A.  Introduction 

This  is  a  report  on  atmospheric  measurements  made  by  the  Naval 
Postgraduate  School  Environmental  Physics  Group  aboard  the  USNS  Kane 
during  operations  in  the  North  Atlantic  in  March  of  1978.  The  primary 
goal  of  these  measurements  was  to  obtain  measurements  of  the  height  of  the 
radar  evaporative  (surface)  duct  from  surface  layer  water  vapor  and 
temperature  profiles  and  to  evaluate  the  practicality  of  using  a  simpli- 
fied (the  so  called  "bulk"  method)  method  for  more  routine  estimations. 
The  bulk  method  requires  measurements  of  only  four  quantities:  the  sea 
surface  temperature  and  the  air  temperature,  wind  velocity  and  relative 
humidity  at  some  reference  height  (usually  10  meters).  It  also  has  the 
advantage  of  considerably  less  stringent  measurement  accuracy  require- 
ments. Due  to  the  extreme  difficulty  in  making  accurate  profile  mea- 
surements at  sea,  the  bulk  method  is  essentially  the  standard  for  evap- 
orative duct  height  calculations. 

B.  Description  of  the  Problem 

In  a  dielectric  medium,  electromagnetic  wave  propagation  is 
influenced  by  the  refractive  index.  Just  as  an  abrupt  change  in  refrac- 
tive index  causes  EM  waves  to  be  "bent",  so  does  a  continuous  gradient  in 
refractive  index.  Since  water  vapor  is  a  major  contributor  to  N,  there 
usually  is  a  strong  gradient  in  N  over  the  sea  surface  because  of  evap- 
oration. This  gradient  causes  EM  waves  to  be  refracted  downward.  If  a 
horizontally  propagating  wave  experiences  a  refraction  greater  than  the 
earth's  curvature,  the  wave  is  said  to  be  ducted.  Since  the  surface  based 
duct  is  caused  by  evaporation,  it  is  often  called  evaporation  duct.  The 
water  vapor  gradient  (and  therefore  the  N  gradient)  over  the  ocean  typi- 
cally decreases  with  height  in  the  surface  layer.  As  a  result,  at  some 
height,  the  N  gradient  has  decreased  below  the  value  required  for  ducting, 
this  is  called  the  evaporation  duct  height,  Z*. 

Since  ducting  can  lead  to  signal  enhancements  of  many  orders  of 
magnitude,  knowledge  of  N  and  Z*  in  the  surface  is  yery   useful.  Fortunately, 
the  physics  of  the  surface  layer  is  well  understood  and  the  specifications 
of  the  N  profile  and  Z*  can  be  reduced  to  the  determination  of  only  four 
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scaling  paramters.  This  theoretical  framework,  known  as  Monin-Obukhov 
Similarity  (MOS),  is  based  upon  the  relation  of  the  surface  layer  profiles 
of  temperature,  wind  speed,  water  vapor  and  turbulence  with  the  surface 
fluxes  of  momentum  and  sensible  and  latent  heat.  As  one  might  expect,  the 
problem  is  simple  in  principle  but  difficult  in  practice,  primarily  because 
the  fluxes  cannot  be  easily  measured  from  an  unstable  platform.  At  present, 
the  best  technique  is  to  measure  the  profiles  at  several  different  heights, 
solve  the  MOS  expressions  "backwards"  to  determine  the  scaling  parameters, 
thus  determining  the  N  profile  and  Z*.  This  process,  equivalent  to  infer- 
ring the  surface  fluxes  from  the  profiles,  is  referred  to  as  profile 
method. 

Multi -level  measurements  necessary  for  the  profile  method  are 
not  trivial.  For  example,  from  measurements  at  heights  of  10  and  25 
meters  above  the  sea  surface,  one  can  typically  expect  differences  on  the 
order  of  10%  in  wind  speed,  3%  in  relative  humidity,  and  0.2C  in  tempera- 
ture. To  determine  the  scaling  parameters  to  a  10%  accuracy,  measurement 
accuracies  required  are:  wind  speed:  1%,  relative  humidity:  0.3%  and 
temperature:  0.02C.  Note  that  the  distortions  caused  by  the  presence  of 
the  ship  could  exceed  these  accuracies.  Atmospheric  measurements  even 
approaching  this  precision  require  highly  trained  scientific  personnel 
with  state-of-the-art  instrumentation.  These  extremely  taxing  measurement 
accuracies  have  led  to  the  interest  in  the  bulk  method. 

If  one  considers  the  MOS  expressions  for  the  gradient  as  mathe- 
matical derivatives,  then  the  integral  of  the  gradient  from  the  surface 
(Z=Z  =0)  to  some  reference  height  (Z)  specifies  the  variable  (say  tem- 
perature, T(Z)).  Note  that  we  have  introduced  two  integration  boundary 
conditions  T(0)  and  Z  ,  where  T(0)  is  the  sea  surface  temperature  and  Z 
is  the  roughness  length.   If  we  assume  that  the  wind  speed  at  the  surface 
is  zero  (U(0)  =  0)  and  the  relative  humidity  at  the  surface  is  100%  (H(0) 
=  100),  only  four  quantities  need  to  be  measured:  T(0),  T(Z),  U(Z),  and 
H(Z).  However,  an  unknown  quantity  Z  ,  which  we  have  also  introduced, 
must  be  empirically  determined. 

Average  temperatures  measured  during  the  cruise  were 
T(0)  =  16. 2C  and  T(Z)  =  14. 7C,  a  difference  of  1.5C  -  an  order  of  mag- 
nitude greater  than  the  two  level  example  given  previously.   In  short,  the 
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bulk  method  has  two  advantages:  fewer  measurements  are  required  and  an 
order  of  magnitude  greater  toleration  to  error.  This  simplicity  was 
obtained  by  introducing  the  quantity  Z  ,  which  has  been  empirically 
studied  by  statistically  averaging  vast  numbers  of  over-water  measurements 
of  bulk,  profile,  and  turbulence  data.  One  of  the  primary  questions  about 
the  bulk  method  is  the  applicability  of  an  "average"  Z  to  a  specific 

situation.  An  even  more  serious  question  deals  with  the  validity  of 

-1 
extrapolating  the  MOS  gradient  expressions  (essentially  Z  )  down  to  the 

surface. 

C.  Summary  of  Measurement  Results 

The  total  data  base  for  the  cruise  consists  of  360  periods  (half 
hour  averages)  of  which  40  periods  contained  near  surface  (Z=2.8m)  bouy 

profile  measurements.  A  tabulation  of  the  EM  propagation  data  (Z*,  Fm, 

2 
N,  dN/dZ  and  C  )  is  given  in  Appendix  C.  There  are  several  ways  of 

looking  at  the  comparison  of  the  profile  or  turbulence  values  with  the 
bulk  values.  For  instance,  the  bouy  profile  values  for  the  water  vapor 
scaling  parameter,  q*,  agreed  wery   well  with  the  bulk  values,  on  the 
average,  but  the  variance  of  individual  values  was  about  50%.  Since  a 
considerable  portion  of  this  50%  is  due  to  uncertainty  in  the  profile 
measurements,  these  results  are  probably  rather  conservative.  The  com- 
parison is  summarized  in  the  following  table: 
QUANTITY       SCALING  PARAMETER  AVERAGE  AGREEMENT 
Water  Vapor         q*  excellent 

Temperature         T*  poor 

Wind  Speed  U*  fair 

Duct  Height         Z*  fair 

D.  Conclusions 
Duct  heights  calculated  on  the  basis  of  bulk  measurements  are 

probably  within  a  factor  of  2  of  the  correct  value,  and  possibly  much 
closer.  The  actual  calculation  of  Z*  from  T(0),  U(Z)  and  H(Z)  can  be 
done  on  a  hand  held  programmable  calculator.  Given  the  sensitivity 
analysis  of  Section  V,  we  feel  that  the  real  problem  in  routine  appli- 
cations is  the  temperature  measurements.  The  wind  speed  can  be  obtained 
with  sufficient  accuracy  using  "eyeball"  estimates  and  a  Beaufort  scale. 


RMS  VARIANCE 

%VARIANCE 

.12  gm/kg 

50% 

.045C 

50% 

.06  m/sec 

15% 

10  M 

100% 
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The  relative  humidity  measurements  can  probably  be  made  with  a  sling 
psychometer  from  a  well  exposed  portion  of  the  ship.  However,  there  is 
no  escaping  the  fact  that  the  air  sea  surface  temperature  difference 
must  be  measured  to  better  than  0.5C.  This  requires  well  calibrated  and 
exposed  temperature  sensors  (with  low  drift  characteristics).  The  air 
temperature  measurement  needs  to  be  made  with  a  good  aspirator  in  a 
location  reasonably  free  of  ship  influence. 

Typical  bulk  Z*  values  estimated  during  the  cruise  were  around 
15  meters.  Under  high  relative  humidity  conditions  (rain  and  fog)  Z* 
was  as  low  as  1  meter.  The  maximum  Z*  was  around  25  meters.  Consid- 
erable time  variations  were  observed,  mostly  associated  with  frontal 
passages.   In  view  of  this,  measurements  should  probably  be  made  several 
times  a  day. 

One  result  of  the  theoretical  analysis  was  the  astounding 
volatility  of  Z*  under  stable  conditions  (air  warmer  than  water).  For 
the  example  with  T(Z)-T(0)=1 .4C,  a  relative  humidity  of  90%  gave  Z*  =  5 
meters  whereas  a  relative  humidity  of  75%  gave  Z*  =   50  meters.  Unstable 
conditions  showed  much  less  variability.   In  fact,  since  the  relative 
humidity  over  the  ocean  is  rarely  less  than  40%,  we  expect  Z*  to  almost 
always  be  less  than  30  meters  during  unstable  conditions.  Therefore, 
stable  conditions  produce  unusually  strong  surface  ducting  and  the  most 
variable  ducting  conditions.  Although  expected,  the  effect  is  larger 
than  we  would  have  anticipated. 

E.   General  Comments  and  Recommendations 

1.  We  feel  another  cruise  is  called  for  to  clear  up  the 
uncertainties  of  the  Kane  data.  We  are  now  "on  line"  with  improved 
temperature  and  humidity  profile  instruments,  an  IR  sea  surface  tem- 
perature sensor  and  an  improved  computerized  data  logger  and  analyzer. 
One  more  cruise  could  do  a  great  deal  in  establishing  the  bulk  method 
procedures. 

2.  An  EM  propagation  analysis  of  other  NPS  data  could  be 
compared  to  the  Kane  data  or  to  climatological  predictions.  We  have  on 
file  data  from  the  Pacific,  the  North  Atlantic  and  the  Mediterranean. 

3.  More  work  needs  to  be  done  on  the  drag  coefficients 
(related  to  Z  ).  The  wind  speed  coefficient  is  fairly  well  known  but 
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real  knowledge  about  the  temperature  and  water  vapor  coefficients  is 
quite  rare. 

4.  Development  of  a  simple  system  for  routine  shipboard 
measurements  of  sufficient  quality  to  allow  bulk  calculations  should 
probably  be  undertaken. 

5.  A  great  deal  of  boundary  layer  work  needs  to  be  done 
(above  the  surface  layer)  to  allow  extension  of  similarity  scaling  to 
greater  heights.   If  successful,  it  may  be  possible  to  make  estimates  of 
elevated  ducts  from  surface  and  remote  sensing  data. 
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II.  Instrumentation 

A.   Sensors  and  Placement 

Rather  than  embark  upon  a  detailed  description  of  the  NPS 
instrumentation,  we  will  include  as  Appendix  A  a  report  (Houlihan  et  al , 
1977)  on  the  subject  and  give  only  a  brief  sketch.  The  mean  measure- 
ments are  summarized  in  the  following  table: 
MEASUREMENT  INSTRUMENT         ACCURACY 

Temperature,  T         Quartz  oscillator         .01 C 
Relative  Humidity,  RH    Li  CI  cell  •  3% 

Wind  Velocity,  V=U      Cup  anemometer  4% 

The  relative  humidity  and  temperature  sensors  are  deployed  in 

aspirated  shelters  to  reduce  errors  from  solar  and  IR  radiation.  These 

aspirators  represent  a  weak  point  in  the  temperature  profile  measurements 

and  do  not  produce  reliable  temperature  profile  data  during  the  day. 

The  turbulence  data  was  based  upon  measurement  of  fluctuations  of  wind 

velocity  (yielding  the  rate  of  dissipation  of  turbulent  kinetic  energy, 

e)  and  temperature  fluctuations  (yielding  the  temperature  structure 

function  parameter,  Cj  ).  Constant  temperature  hot  film  sensors  were 

used  to  measure  e   and  microthermal  platinum  wires  were  used  to  measure 

2 
Cj  .  The  frequency  response  of  these  measurements  is  greater  than 

500HZ. 

The  instruments  are  deployed  at  discrete  levels:  two  levels 

on  the  NPS  tower,  one  on  the  bow  davit  and  (circumstances  permitting) 

one  on  the  buoy  (see  Figure  la).  The  following  are  the  level  heights 

(Z)  above  the  water: 

Level  #  12        3        4 

Height,  meters     2.8      9.1      18.6     24.7 
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Z1  =  2.8m 
Z2=  9.1m 
Z3=  18.6m 
Z4=24.7m 


Sea  Surface  Temperature 


BOUY 


Figure  la.  Mounting  Arrangement  on  USNS  Kane  for  NPS  Equipment. 

The  bouy,  which  is  lowered  into  the  water  with  a  crane  on  the 
bow  of  the  ship,  is  used  to  obtain  a  measurement  very   close  to  the 
surface.  (See  Figure  lb).  Since  the  gradients  we  wish  to  measure  vary 
roughly  as  Z"  ,  they  are  easier  to  measure  near  the  surface.  Unfortu- 
nately, the  buoy  can  only  be  used  when  a)  the  ship  is  stationary, 
b)  the  seas  are  moderate. 
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Figure  lb.  The  Bouy  Deployed  from  USNS  Kane 
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III.  Theory 

A.   EM  Propagation 

1.  Refractivity  N  and  Ducting 

The  refractivity  N  (in  N  units)  of  the  atmosphere  for  EM 
wave  propagation  is 

N  =  77.6  P/T  +  3.73  x  105  e/T2  (1) 

Where  P  is  the  atmospheric  pressure,  T  the  absolute  temperature  and  e  is 
the  partial  pressure  of  water  vapor.  The  vertical  gradient  of  N  is 

dN  =  c  IE.  +  c  M+  c  II  (?) 

dZ   Ll  3Z   L2  3Z   L3  3Z  {Z) 

where  q  is  the  water  vapor  mixing  ratio  (q  =  .625  e)  in  gm/kg  and 
C,  =  .3,  C~  -  7.2  and  C3  -  -1.3.  A  negative  gradient  in  refractivity  at 
height  7.  causes  a  horizontally  propagating  EM  wave  to  be  bent  towards  the 
surface.  If  the  gradient  is  large  enough,  the  amount  of  bending  will  ex- 
ceen  the  curvature  of  the  earth  and  the  wave  will  be  ducted  at  the  height 
Z.  Since  the  magnitude  of  the  gradient  decreases  with  increasing  height 
in  the  surface  layer,  at  some  height  (Z*)  the  gradient  will  no  longer  be 
strong  enough  for  ducting.  Waves  propagating  below  Z*  are  ducted,  those 
above  are  not.  This  height  is  called  the  evaporative  duct  height,  and 
is  defined  by  the  critical  gradient  necessary  for  EM  wave  trapping. 

dN  _   ,C7  -1  ,-N 

dl  "  --157  m  (3) 

Using  Eq.  21  and  noting  that  dP/dZ  =  -pg  (the  normal  hydrostatic  balance 
for  the  atmosphere),  the  Z*  is  defined  by  the  following  equation 

-.157  =  -.032  +  C2if+C3|I  (4) 

Since  3q/3Z  and  3T/3Z  are  height  dependent  (approximately  ^Z"  )  there 
will  exist  a  value  of  Z  =  Z*  such  that  Eq.  4  is  satisfied. 
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2 

2.  Structure  Function  C., 

Turbulence  in  the  atmosphere  produces  fluctuations  in  N  on 
a  broad  range  of  size  scales  (milimeters  to  kilometers).  These  fluctua- 
tions scatter  the  EM  waves  producing  a  loss  of  resolution  and  coherence. 
If  N(x)  and  N(x+d)  represent  values  of  N  at  two  points,  separated  a  dis- 
tance d,  then 

CN2  =  <(N(x)  -  N(x+d))2>d"2/3  (5) 

For  EM  Waves 

CN2  ■  <  3  '  2  ^  +  I  (  '  +  1  T  »  CTq3  (6) 

2 

where  A  and  B  are  constants  C   is  the  water  vapor  structure  function 

q 

and  C-j-  is  the  covariance  function.  Neglecting  the  Cy  term  (at  sea  level) 
C  2  =  [  (77.6  x  10"6)(7.67  x  103)-,2  ,y) 

B.   Monin-Obukhov  Similarity  Surface  Layer  Equations 

The  vertical  gradients  of  the  wind  speed  (U),  potential  tem- 
perature (T)  and  water  vapor  mixing  ratio  (q)  are  given  by 

dX/dZ  ■  ^Z     *x  <«)  <8' 

Where  X  =  U,  T  or  q,  <   =   Von  Karmon's  constant  (35),  a  is  the  diffusi- 

A 

vity  constant  (a  =1.0),  X*  is  the  scaling  parameter  and  <j>(c)  is  the 
MOS  stability  correction  function  (Appendix  B).  Also 

I   =  Z/L  (9) 

where  L  is  the  Monin-Obukhov  stability  length  defined  as 

L  =  .9(t/"*J8  q.)  HO) 

If  we  integrate  Eq.  8  from  Z  =  lnv   to  Z,  we  obtain 

ox 

X(Z>  =  X<Zox»+rt    nnZ/Zox-*x  (5)]  (11) 

A 

Where  iJ>Y(0   is  the  profile  stability  function   (Appendix  B). 

2 
The  structure  function  parameter  C       is  related  to  X*  by  the 

A 

function  fx(c)  (Appendix  B) 

Cx2  =  X*2  Z'2/3  fx  U)  (12) 

A  similar  form  exists  for  e  ,  the  rate  of  dissipation  of  fluctuations 
of  X 

S  =  BxCx2  *y2  =  L^tW  (") 

Where  6  is  the  Corrsin  constant  and  E  (?)  is  the  MOS  dissipation  func- 

A  A 

tion.   In  this  format,  e  =  e. 
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C.   MOS  Parameter  Determination  Methods 
1 .   Profile  Method 

The  profile  method  is  based  upon  Eq3  11,  expressed  in  the 
following  form 

x(z)  =^[m  z-*x(0]  ♦  [x(0  -  ^  m  zox]       (i4) 

A  X 

This  can  be  expressed  as  a  linear  regression  of  the  form 

X(Z)  ■  mY  +  b  (15) 

Where  Y  =  In  Z  \\>   (s).  If  we  have  two  or  more  levels  of  X(Z),  we  can  fit 
the  date  to  Eq.  16,  the  slope  m,  of  the  fit  gives  X*  =  a  <m.  Note  that 

A 

it  is  not  necessary  to  know  the  constant  factor,  b,  (related  to  Z  ) 
to  find  X*.  The  actual  process  is  done  iteratively  in  the  following 
manner: 

a.  Assume  a  value  for  L 

b.  For  each  data  point  at  height  Z,  calculate  £  and 

♦x(e). 

c.  Least  squares  fit  the  multi-level  data  to  obtain 
U*,  T*  and  q*. 

d.  Calculate  a  new  value  of  L  from  Eq.  10 

e.  Return  to  b)  until  L  converges 

Since  the  wind  speed  profiles  from  the  Kane  were  consi- 
derably influenced  by  the  ship,  we  used  the  bulk  method  to  calculate 

u*. 

2o   Bulk  Method 

The  bulk  method  is  also  based  on  Eq.  15.  If  one  knows 
values  for  ZnY,  then  X*  can  be  calculated  simply  from  X(Z)  and  X(0). 

UA 

Rewriting  Eq.  15 

y  -      (X(Z)  -  X(0))  nc. 

X*  "  V  (In  Z/ZQx  - Xti))  <16) 

We  see  the  form  of  the  standard  drag  coefficient  equation 

X*  =  Cx1/2  (X(Z)  -  X(0))  (17) 


Where  C»  is  the  drag  coefficient  for  X 


r 

XN 

X    ~  7t  Tin  o"  (18) 

[i  -  (V)  ]  cXN1/2  ♦(5)]2 
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In  this  form,  Z   is  related  to  the  neutral  drag  cofficient  by 

1/2      axK 
CXN    =  In  (Z/ZQx)  (19) 

The  MOS  stability  parameter,  £,  can  be  obtained  from 
Eq.  10  and  Eq.  18. 

,  =  ,Q     []  l£]  CUN1/2^u^))2  (20) 

(1  -  (a^)"1  CTN1/2  ^(?)) 

Where  -..p 

5  s<gZ  ^TN [(T(Z)  -  T(0))  +  .18(g(Z)  -  g(0))]     (21) 

0    T    CUN  U(Z)2 

In  this  case,  the  process  is  somewhat  simpler,  one  cal- 
culates e  from  the  data  and  solves  Eq.  21  iteratively  for  £.  Given  £, 
the  MOS  parameters  are  calculated  straightforwardly  from  Eq.s  18  and  19. 
3.   Turbulence  Method 

Given  a  reasonable  estimate  of  %   (from  either  of  the  two 
previous  methods)  one  can  obtain  an  independent  value  of  the  scaling 
paramete 
(Eq.  12) 


2    2 
parameter  from  turbulence  measurements  of  CT  ,  C   or  e.  For  instance 

T  q 


T*2  =  CT2  Z2/3/fT(?)  (22) 

or   (Eq.    13) 

U/  =  eZ/EjjCc)  (23) 

D.   EM  Propagation  Equations  in  MOS  Form 

Near  sea  level,  the  gradient  of  M  can  be  written  (Eq.s  2  and  3) 

*TU) 

dN/dZ  =  -.032  +  ^y-  (7.2  q*  -  1.3  Tj  (24) 

The  duct  height,  Z*,  is  obtained  from  Eq.  25  by  setting 
dN/dZ  =  -.157,  the  critical  gradient  for  ducting. 

-(7.2  q*  -  1.3  TJ 
Z*  =  " aT<  (.125)    ¥Z*/D  (25) 

Given  q*,  T*  and  L,  we  solve  Eq.  26  iteratively  to  obtain  Z*.  Note 
that  no  ducting  will  occur  unless 

7.2  q*  -  1.3  T*  <0  (26) 
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The  refractive  index  structure  parameter  is  calculated  from  Eq.s  7 
and  12 

CN2S[(77-6X10"6)(7'67X1°3)]2  q*2  Z"2/3  fT(C)      (27) 

IV.  RESULTS 

A.   Discussion 

Before  presenting  the  results,  it  is  worthwhile  to  consider 
several  points  about  the  absolute  validity  of  the  entire  process.  The 
various  stability  correction  functions  -  $(€),  f(0»  ^(0  -  are  empiri- 
cally determined  functions.  In  fact,  they  are  based  on  overland  mea- 
surements. The  NPS  group  has  spent  a  great  deal  of  effort  in  evaluating 
the  validity  of  these  functions  over  water  and,  so  far,  they  appear 
to  be  quite  good.  There  are  other  problems  with  some  of  the  empirical 
constants  and  the  drag  coefficients.  The  profile  method  is  dependent 
upon  the  product  cu-ic,  both  measured  empirically.  The  estimates  we 
have  used  are  a_  =  1.35  and  <   =  .35,  from  Businger  et  al  (1971).  The 
actual  value  of  both  constants  is  still  a  subject  of  controversy. 

The  bulk  method  is  dependent  upon  the  drag  coefficients 
C,,N,  CTN  and  C  ...  For  our  analysis,  we  assumed  that  water  vapor  has 
the  same  characteristics  as  temperature,  that  is 

aT  =  aq  *'  Z0T  =  Z0q  ;  *T(c)  =  *q(0 
There  are  theoretical  reasons  to  believe  that  this  is  a  good  approxi- 
mation but  not  exactly  correct.  Under  these  assumptions,  CjN  =  C  ... 
For  the  actual  form  of  the  drag  coefficients,  we  assumed  C.JN  had  the 

wind  speed  dependence  given  by  J.  Kondo  (1975).  Based  upon  measure- 

2 
ments  summarized  by  Liu  (1977)  and  NPS  measurements  of  CT  ,  we  assigned 

-3 
a  value  of  CTN  =  1.1  x  10   and  solved  for  c  as  a  function  of  5  . 

Thus,  for  each  data  value  of  5  we  directly  calculated  5  from  the  for- 
mula 

5  =  50(1  +  .IV4)   ?Q<0  (28a) 

5  =C0d  +  .^0-4+  .26CQ2J)  ?o>0  (28b) 
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B.   EM  Propagation  Results 

A  listing  of  the  bulk  calculations  of  the  EM  data  is  given  in 
Appendix  C  (Table  I-C).  The  data  includes  calculations  of  the  Z=10 
meter  values  of  N,  dN/dZ  and  CN  as  well  as  Z*  and  Fm.  Fm  is  the  mini- 
mum frequency  trapped  by  the  duct  as  is  calculated  from 


'  3.6  x  102  Z*"3/2 


(29) 


Where  Fm  is  in  GHz. 

The  following  are  two  examples  to  illustrate  the  general  be- 
havior of  Z*.  Figure  2  shows  Z*  for  a  three  day  period  at  anchor  ap- 
proximately 80  miles  west  of  the  Straits  of  Gibraltar.  The  steady  de- 
crease of  Z*  for  this  first  part  of  the  record  was  associated  with  high 
relative  humidity  and  rain.  The  rapid  increase  of  Z*  around  0200  hours 
on  3/3  corresponded  to  the  passage  of  the  front  as  indicated  by  clear- 
ing skies  and  a  shift  of  the  wind  from  210°  to  260°.  Figure  3  shows 
a  very   similar  event  on  station  in  the  survey  area.  In  this  case,  the 
weather  system  was  not  strong  enough  to  produce  rain  but  fog  was  ob- 
served at  1800  hours  on  3/10.  The  wind  shifted  from  200°  to  300  be- 
tween 1900  hours  on  3/10  and  0700  hours  on  3/11. 
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Figure  2.  Z*  values  for  period  3/1  to  3/4  when  KANE  was  anchored  80  miles 
west  of  Straits  of  Gibraltar. 
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Figure  3.  Z*  values  for  period  3/10  to  3/12  when  KANE  was  in  survey  area. 

C.   Evaluation  of  the  Bulk  Method 
1 .  Seal ing  Parameters 

In  order  to  make  this  evaluation  more  meaningful,  we  have 
analyzed  the  scaling  parameters  under  those  conditions  of  maximum  profile 
accuracy.  For  instance,  the  humidity  measurement  were  not  accurate  enough 
to  give  reliable  profiles  except  when  the  bouy  was  deployed.  The  tempera- 
ture profiles  were  only  analyzed  for  the  night  time  periods.  The  turbu- 
lence measurements  of  U*  were  valid  during  all  periods.  A  complete  tabu- 
lation of  comparisons  for  all  periods  is  given  in  Table  II-C  (Appendix  C) 
but  keep  in  mind  that  the  graphs  shown  below  are  based  on  data  selected 
from  the  table.  Incidently,  the  bouy  data  comparison  has  been  extracted 
for  condensed  viewing  and  is  shown  in  Table  I II-C  (Appendix  C). 

The  comparison  of  bulk  and  profile  calculations  of  q*  are 
shown  in  Figure  4  from  the  bouy  data.  On  average,  the  comparison  is 
yery   good.  The  individual  point  scatter  is  about  50%.  For  this  graph, 
and  those  following,  the  error  bars  are  for  the  mean  estimates  and  the 
number  represents  the  number  of  data  points  in  each  bin. 
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Figure  4.  Calculation  of  water  vapor  mixing  ratio  scaling  parameter, 
q*,  profile  determination  vs.  bulk,  bouy  data  only. 

The  comparison  of  bulk  and  profile  calculations  of  T*  are  shown 
in  Figure  5.  These  results  are  very   poor  indeed.  It  is  not  clear  if 
the  lack  of  correlation  is  due  to  measurement  problems  or  is  in  fact 
failure  of  the  bulk  method. 
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Figure  5.  Calculation  of  potential  temperature  scaling  parameter,  T*, 
profile  determination  vs.  bulk,  night  time  data. 

The  comparison  of  the  bulk  and  turbulence  calculations  of  U*  are 
shown  in  Figure  6.  Although  the  correlation  is  not  perfect,  it  is  pro- 
bably good  enough  for  bulk  determinations  of  Z*.  The  complicated  velo- 
city dependence  of  CUN  that  was  used  for  these  U*  bulk  values  did  not 

give  significantly  better  comparison  than  assuming  a  velocity  indepen- 

_3 
dent  value  of  C.]N  =  1.3  x  10  . 
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Figure  6.  Calculation  of  friction  velocity,  U*,  turbulence  determina- 
tion vs.  bulk,  all  data. 


2.   Duct  Height 

Evaporation  duct  height,  Z*,  was  calculated  using  the 
bulk  and  profile  methods  for  all  data.  A  comparison  of  these  results 
is  given  in  Figure  7. 


28 


1.0 


2.0 


4.0  6.0   10    20 
BULK  Z«  (m) 


40  60   100 


Figure  7.  Calculation  of  evaporation  duct  height,  Z  ,  profile  deter- 

A 

mination  vs.  bulk.  The  open  circles  represent  all  data,  the  circled 
x's  are  the  bouy  data  only. 

Although  the  correlation  is  only  fair,  it  should  be  noted 
that  the  bins  containing  most  of  the  points  (187  out  of  232)  fall 
within  the  expected  mean  error.  Also  encouraging  is  the  much  better 
agreement  obtained  from  the  bouy  data. 
V.   SENSITIVITY  ANALYSIS 
A.   Data  Distribution 

Rather  than  consider  many  possible  combinations  of  the  me- 
teorological parameters,  the  sensitivity  analysis  will  be  done  with- 
in the  framework  of  typical  values  measured  during  the  cruise.  A 
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summary  of  the  overall  average  data  is  given  below. 

QUANTITY  SYMBOL      AVERAGE 

Sea  Surface  Temp.       TS        16.25°  C 
Air  Temperature        T(10)      14.67°  C 
Relative  Humidity       H(10)      70% 
Wind  Speed  U(10)      7.57  m/sec 

Mixing  Ratio  q(10)      7.29  gm/kg 

B.   Sensitivity  of  Z* 

The  sensitivity  analysis  was  performed  using  the  meteorologi- 
cal parameters  as  inputs  into  the  bulk  expressions.  The  predictions 
are  shown  in  the  form  of  Z*  as  a  function  of  relative  humidity.  Of 
course,  one  should  keep  in  mind  the  extreme  rarity  of  relative  humidi- 
ties (at  Z=10  meters)  less  than  40%  over  the  ocean.  On  each  of  the 
graphs  shown,  the  solid  line  represents  the  "average"  curve  (TS-16.2, 
T(10)-14.7  and  L)(10)=7.6  m/sec).  Since  the  average  relative  humidity 
was  70%,  we  expect  an  average  Z*  =  16  meters.  Figure  8  shows  the  ef- 
fects of  velocity  variations  on  the  order  of  the  standard  deviation 
found  during  the  cruise.  The  effect  of  sea  surface  variation  (at  con- 
stant air-sea  temperature  difference)  is  given  in  Figure  8b.  Z*  is 
moderately  insentive  to  both  parameters.  The  corresponding  graph  for 
sensitivity  to  the  air-sea  temperature  difference  is  noticeably  differ- 
ent (Figure  9).  Not  only  can  \/ery   large  values  of  Z*  occur  for  stable 
conditions  (AT>0)  but  note  the  yery   steep  variability  of  Z*  between 
60%  and  90%  humidity.  Of  course,  one  must  bear  in  mind  that  these  are 
surface  layer  extrapolations,  the  maximum  height  of  validity  is  proba- 
bly less  than  50  meters.  Nevertheless,  these  results  are  quite  sig- 
nificant. 


30 


100 
70 

50 

40 

30 
20 


N 


10 
7h 

5 

4 

3  - 
2  - 


_ 

- 

Ts= 

16.2  °C 

- 

AT= 

-1.5  ■C 

- 

^ 

U=10.4  m/sec 

- 

u= 

=4.8  m/sec      ^-x    ^^ 
'n 

7.6  m/sec 
NX 

- 

- 

\ 

- 

__1_ 

III! 

l 

!             I             I 

0     10  20  30   40  50  60   70   80   90  100 
RH(%)  8a 


100 
70 
50 

30 
20 


% 

N 


10 


5  - 
4- 


T-V-!5°C 
U=7.6  m/sec 


T>18.2°C 


TS=16.2°C 


T5=14.2°C 


i       i       i       '       '       i       i       i 


0     10    20  30  40  50  60   70   80  90  100 
RH(%)  8b 


Figure  8.     Duct  height  vs.   relative  humidity  for  typical  meteorological 
conditions  found  during  the  Kane  cruise;   a)   sensitivity  to  wind  speed 
b)   sensitivity  to  sea  surface  temperature. 
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Figure  9.     Duct  height  vs.  relative  humidity  for  typical  meteorological 
conditions  found  during  the  Kane  cruise:     sensitivity  to  the  air-sea 
temperature  difference  .(AT  =  T(Z)  -  T(0)). 
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The  large  difference  in  stable  and  unstable  atmospheric 
evaporation  ducting  properties  is  a  direct  result  of  the  behavior  of 
the  MOS  gradient  equations.  Let  us  assume  that  yery   near  the  surface 
we  have  an  N  gradient  great  enough  to  permit  ducting.  Since  the  gra- 
dient very   near  the  surfaces  decreases  with  height  as  Z  ,  we  expect 
the  gradient  to  decrease  with  increasing  height  until  Z=Z*  and  the  cri- 
tical gradient  is  reached.  The  atmospheric  stability  enters  into  the 
problem  through  its  influence  on  the  height  dependence  of  the  gradients. 

Under  unstable  conditions,  the  convective  mixing  quickly  removes  the 

-3/2 
gradients  as  Z  increases,  leading  to  a  Z    gradient  height  dependence. 

Consequently,  the  critical  gradient  is  reached  at  a  relatively  low  Z. 
Under  stable  conditions,  the  bouyancy  forces  are  retarding  the  turbu- 
lent mixing  as  Z  increases,  leading  to  gradients  that  approach  a  con- 
stant independent  of  height.  Therefore,  the  gradient  decreases  more 
slowly  under  stable  conditions  and  the  critical  gradient  can  be  reached 
at  much  greater  heights. 
VI.  CONCLUSIONS 

A.  Variability  of  Z* 

Typical  values  of  Z*  obtained  during  the  cruise  ranged  from 
around  2.0  meters  to  25  meters  with  a  mean  value  of  15  meters.  Trans- 
missions at  10  GHz  would  have  been  ducted  approximately  70%  of  the 
time,  20  GHz  approximately  85%  of  the  time.  The  two  factors  having 
the  greatest  influence  on  Z*  are  the  relative  humidity  and  the  air-sea 
temperature  difference.  In  one  instance  the  duct  height  changed  from 
3.0  meters  to  20  meters  in  about  four  hours  following  a  frontal  passage. 

B.  Validity  of  the  Bulk  Method 

The  profile  and  turbulence  data  from  this  cruise  do  not  give 
a  clear  validation  of  the  bulk  method.  Although  the  water  vapor  and 
wind  speed  scaling  parameters  were  reasonably  well  correlated  in  the 
mean,  there  was  considerable  scatter  (on  the  order  of  50%  and  15%  res- 
pectively) for  individual  measurements.  Comparisons  for  T*  and  Z*  were 
definitely  a  disappointment.  We  feel  that  much  of  the  discrepancy  can 
be  attributed  to  inaccurate  temperature  and  humidity  profiles.  Since 
both  of  these  systems  have  been  upgraded  winch  the  Kane  cruise,  another 
cruise  should  be  made.  Based  upon  comparison  of  turbulence  T*  (from  CT  ) 
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with  bulk  T*  from  other  NPS  data,  we  feel  that  the  bulk  T*  are  much  better 
than  indicated  by  the  Kane  data. 

C.   Implications  of  the  Sensitivity  Analysis 

The  sensitivity  analysis  showed  that  the  wind  speed  and  sea  sur- 
face temperature  contributions  to  Z*,  though  significant,  were  not  highly 
critical.  The  relative  humidity  and  air-sea  temperature  difference  were 
definitely  critical  factors.  Given  the  standard  measurement  capabilities 
of  these  two  quantities,  it  is  obvious  that  the  air-sea  temperature  dif- 
ference will  be  the  primary  roadblock  to  implementing  a  program  of  rou- 
tine measurement  of  Z*  either  from  ship  data  or  satellite  data.  Another 
result  of  the  study  was  the  large  values  of  Z*  possible  under  stable  con- 
ditions (stable  conditions  predominate  off  the  east  coast  of  the  U.S.  and 
Canada  in  the  summer). 
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ABSTRACT 

Programs  of  investigation  of  optical  propagation,  marine  fog,  and 
boundary  layer  properties  in  the  marine  environment  have  led  to  the 
development  of  a  shipboard  system  for  measuring  atmospheric  stability 
and  turbulence.   This  research  requires  very  accurate  multilevel 
measurements  of  mean  and  fluctuating  temperature,  wind  velocity,  and 
humidity.   This  paper  is  a  description  of  the  system  including  the 
following  aspects:   1)  instrumentation,  2)  calibration,  3)  analysis 
techniques,  and  4)  examination  of  the  data  and  evaluation  of  the 
system. 
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I.   Introduction 

There  is  currently  considerable  interest  in  turoulent  transport  pro- 
cesses in  the  marine  boundary  layer.   Two  problems  which  are  of  parti- 
cular importance  concern  some  turbulence  properties  and  their  influence 
on  optical  propagation,  and  the  correlation  between  the  transport  of  heat 
and  moisture  and  the  occurrence  of  marine  fog.   In  the  optical  problem, 
turbulence  causes  beam  wander  and  spreading  and  scintillation  (twinkling) , 
which  cause  energy  density  loss,  resulting  in  information  loss  in  the  case 
of  communication  systems.   The  occurrence  of  marine  fog  is  strongly  depen- 
dent on  the  moisture/heat  ratio  in  the  atmosphere  and  the  turbulent  trans- 
port of  these  quantities  is  of  critical  importance.   The  goal  of  present 
investigations  of  turbulence  properties  is  to  enable  predictions  of  the 
relevant  atmospheric  properties  on  the  basis  of  relatively  simple  obser- 
vations that  are  currently  being  used  by  weather  centers.   A  sophisticated 
series  of  experiments  is  necessary  at  this  stage  of  development  because 
of  the  need  to  develop  an  adequate  parameterization  of  the  problem  and 
because  of  the  complications  inherent  in  model  development  and  verification. 

The  Naval  Postgraduate  School  research  vessel  R/V  Acania  is  the 
primary  platform  for  overwater  measurements  and  it  is  the  system  which  is 
installed  on  this  ship  which  is  described  here.   A  four-level  system  has 
been  developed  to  measure  both  mean  and  fluctuating  quantities.   This 
system  is  augmented  by  a  visiometer,  a  sea  surface  temperature  probe,  an 
acoustic  sounder,  and  at  times  by  balloons,  kytoons  or  kites  to  make  ele- 
vated measurements. 

Since  a  ship  is  a  non-stable  platform  there  is  always  some  question 
as  to  the  validity  of  data  obtained  when  turbulence  parameters  are  being 
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investigated.   It  is  not  possible  to  measure  turbulence  directly  at  low 
frequencies  (<lHz)  aboard  ship  because  the  ship  motion  interferes  in 
this  region.   At  higher  frequencies  valid  data  can  be  obtained,  and  the 
experimental  and  analysis  techniques  described  herein  exploit  this  fact. 
An  additional  problem  with  shipboard  measurements  is  the  large  physical 
size  of  a  ship  and  the  disturbance  of  the  local  air  flow  by  the  ship's 
presence.   Fortunately,  the  R/V  Acania  is  a  fairly  small  and  narrow  ship 
and  this  effect  is  not  large.   Moreover,  proper  placement  of  sensors  to 
take  data  only  when  the  wind  direction  is  favorable  can  minimize  inter- 
ference effects  on  the  data.   In  this  paper  evidence  is  presented  to  show 
that  neither  ship  motion  nor  physical  presence  invalidate  the  data 
displayed . 


II .   Shipboard  Configuration 

The  shipboard  sensor  mounting  arrangement  is  shown  in  Figure  1.   The 
sensors  for  mean  and  fluctuating  wind,  temperature  and  humidity  were 
mounted  on  two  towers,  located  approximately  2  feet  and  16  feet  behind 
the  tip  of  the  bow.   The  heights  of  the  sensors  above  the  sea  surface  were: 

Level  1  4.2  meters 

Level  2  6.6  meters 

Level  3  7.6  meters 

Level  4  13.9  meters 

The  following  sensors  were  installed  at  the  levels  indicated  on  the  masts: 


Mean  Wind 

Mean  Temperature 

Mean  Humidity 

Wind  Fluctuation 

Temperature  Fluctuation 

Humidity  Fluctuation 


Cup  Anemometer 
Quartz  Thermometer 
Li  CI  Cell 
Hot  Wire 
Cold  Wire 
Lyman  -  a 


All  Levels 
All  Levels 
All  Levels 
All  Levels 
All  Levels 
Levels  1  and  3 
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The  heights  and  locations  of  the  sensors  were  chosen  in  such  a  way 
to  minimize  shipboard  influence  and  to  measure  over  a  wide  enough  height 
difference  so  that  small  gradients  in  the  parameters  would  be  detected. 
A  subsequent  improvement  that  allows  data  to  be  obtained  much  closer  to 
the  sea  surface  will  be  described  at  the  end  of  this  section. 

The  acoustic  sounder  was  mounted  on  the  fantail  of  the  ship  since 
this  is  the  area  that  is  most  noise  free.   The  vibrations  of  the  deck 
when  the  ship  is  underway  produce  enough  noise  that  useful  signals  can 
only  be  obtained  up  to  about  400  meters. 

The  visiometer  was  mounted  on  top  of  the  bridge  on  a  3  feet  high 
support.   This  position  placed  the  instrument  in  an  uncontaminated  air- 
flow when  the  relative  wind  was  within  90  of  the  bow. 

The  sea  surface  temperature  sensor  was  located  approximately  4  feet 
from  the  starboard  side  of  the  ship  and  floats  within  one  inch  of  the 
surface.   The  starboard  side  was  chosen  because  the  ship  water  exhausts 
were  located  on  the  port  side. 

In  order  to  study  near  surface  effects  and  to  take  advantage  of 
larger  gradients  in  the  near  surface  region,  a  buoy  system  (Figure  2) 
was  developed  for  deployment  off  the  bow  of  the  Acania  (Corbin,  1977) . 
A  heavy  weight  was  suspended  by  cable  through  a  hollow  mast  in  the  buoy, 
thus  allowing  the  buoy  to  be  rigidly  controlled  in  the  horizontal  dir- 
ection but  to  freely  slide  up  and  down  with  waves. 

III.   Instrumentat  ion 
A.   Mean  System 

A  block  diagram  of  the  mean  system  is  shown  in  Figure  3.  The 
system  consisted  of  two  primary  components,  viz.,  the  sensors  and  their 
associated  electronics,  and  the  data  logging  and  readout. 
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The  data  logger  was  a  microprogrammable  integrated  data  acquisition 
system  (MIDAS)  which  was  developed  at  NPS  specifically  for  acquiring 
shipboard  meteorological  data.   MIDAS  utilized  an  Intel  Model  8008 
Central  Processor  to  control  the  sampling,  averaging  and  recording  of 
mean  meteorological  data.   All  software  programs  were  written  in  PL/M 
to  facilitate  the  writing  of  self -documenting  programs  (Plunkett,  1976) . 

The  operator  was  interfaced  with  the  system  via  teletype  for  full 
duplex  input/output  communications  and  program  control.   The  operator 
exercised  control  over  the  sample  start-time  and  the  number  of  samples 
to  be  averaged  before  outputting.   The  operator  could  also  alter  the 
preset  sample  list  by  adding  or  deleting  various  sensors  as  they  came 
on  line  or  become  inoperable.   Once  initiated,  the  system  was  fully 
automated  to  sample  the  tailored  list  of  sensors  every  30  seconds  and 
periodically  printed  output  values  averaged  over  a  selected  interval 
from  one  to  sixty  minutes. 

Output  values  were  printed  on  the  teletype  in  columnized  format 
with  the  time  of  print  as  a  leader.   The  teletype  contained  a  paper 
tape  punch  that  could  be  activated  by  the  operator  to  produce  a  data 
copy  concurrent  with  the  printout.   A  magnetic  cassette  tape  recorder 
was  integrated  into  the  system  as  a  third  data  output  device.   This 
cassette  is  to  be  interfaced  to  a  Hewlett-Packard  Model  9831  Portable 
Computer  so  that  profile  and  gradient  flux  estimates  can  be  performed 
onboard  automatically,  using  the  BASIC  programming  capability  of  the 
HP  9831.   Presently  data  cards  are  punched  from  the  paper  tape  output 
and  processed  on  the  IBM  360  System  at  NPS. 

Mean  wind  speeds  were  measured  using  Thornthwaite  Cup  Anemometer 
units  featuring  lightweight  plastic  sensors  that  ensured  low  wind  speed 
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response.   The  pulse  outputs  (one  pulse  per  revolution)  of  the  four 
anemometer  units  were  fed  directly  to  MIDAS  where  they  were  averaged 
and  printed  on  the  teletype  at  the  conclusion  of  each  preselected  time 
period.   Revolution  counts  were  also  displayed  on  Hewlett-Packard  Model 
5221A  Counters  for  quick-look  and  calibration  checks.   The  counts  were 
converted  to  wind  speed  using  the  calibration  parameters  supplied  by 
Thornthwaite,  which  were  checked  periodically  in  a  wind  tunnel. 

The  temperature  trnasducers  utilized  were  Hewlett-Packard  Model  2833 
Oceanographic  (quartz  crystal)  Sensors,  with  a  resolution  of  0.01  C. 
The  sensors  and  accompanying  booster  oscillators  were  sequentially  multi- 
plexed through  a  Hewlett-Packard  Model  2801  Quartz  Thermometer  Readout 
Unit  for  monitoring. 

The  humidity  sensors  utilized  were  Hygrodynamics,  Inc.,  Model  1818W 
Wide  Range  Sensors,  which  featured  an  integral  thermistor  probe  with  a 
resolution  of  0.1  C.   These  units  were  sequentially  multiplexed  through 
a  Hygrodynamics  Digital  II  Readout  Unit  for  monitoring  both  relative 
humidity  and  temperature  data. 

The  temperature  probes  in  the  humidity  sensor  arrangements  provided 
a  check  capability  on  the  more  sensitive  quartz  thermometer  sensors 
colocated  with  the  humidity  sensors.   While  the  resolution  demanded  for 
gradient  measurements  (0.01  C)  could  only  be  satisfied  by  the  quartz 
thermometer  system,  the  availability  of  back-up,  bulk  temperatures  from 
the  humidity  units  was  useful  on  several  occasions  where  conflicting 
profile  trends  were  discerned  in  quartz  temperature  readings. 

Temperature  and  humidity  sensor  outputs  at  each  of  the  four 
measurement  levels  were  multiplexed  to  their  respective  readout  units 
using  an  NPS  developed  level  selector  which  operated  under  the  timing 
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control  of  the  MIDAS  system.   The  level  selector  fulfilled  two  essential 
requirements.   First,  it  allowed  one  readout  unit  to  serve  the  entire 
profile  array,  thus  decreasing  the  capital  investment  in  equipment. 
Secondly,  the  level  select  unit  contained  the  power  elements  for  each 
level  in  the  array.   Thus,  when  individual  sensors  were  selectively 
connected  to  their  respective  readout  unit,  only  signal  levels  were 
switched  at  each  junction.   Maintaining  constant  powering  at  each  level 
rather  than  switching  in  power  concomitant  with  signal  level  switching 
kept  line  transients  to  a  minimum  and  insured  proper  time  constant 
response  for  the  temperature  and  humidity  equipment. 

The  quartz  temperature  sensors  and  oscillator  circuit,  and  the 
humidity  sensors  at  each  level  were  contained  within  C.C.  BREIDERT  CO., 
Air-X-Houster  Aspirated  Weather  Shelters.   The  mounting  arrangement  is 
shown  in  Figure  4.   The  shelters  served  two  purposes,  viz.,  protecting 
the  sensors  from  the  marine  environment,  and  reducing  radiation  effects. 
Radiation  effects  are  especially  severe  in  this  application  since  some 
of  the  sensors  are  above  the  water  and  some  above  the  deck,  of  the  ship, 
which  is  heated  by  the  sun  and  radiates  strongly.   Because  of  radiation 
from  the  deck  the  shelters  were  modified  by  placing  a  shield  at  the  base 
of  the  unit  to  reduce  radiation  from  below. 

B.   Fluctuation  System 

A  block  diagram  of  the  fluctuation  system  is  shown  in  Figure  5. 
It  is  immediately  apparent  that  there  are  several  output  devices,  allow- 
ing several  methods  of  data  analysis.   The  analog  tape  recorder  was  used  to 
record  the  basic  signals,  which  could  then  be  analyzed  in  the  laboratory 
subsequent  to  a  field  trip.   The  analog-digital  converters  and  the 
averaging  and  storage  were  contained  in  the  MIDAS  system  described  above 
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and  their  results  were  printed  on  the  teletype  at  the  same  time  as  the 
mean  system  data.   Averaging  could  also  be  accomplished  directly  on  the 
strip  chart  records. 

All  sensors  were  mounted  on  wind  vanes,  with  one  vane  per  level, 
to  maintain  correct  alignment  with  the  relative  wind  directions. 

The  sensors  for  wind  speed  fluctuations  were  Thermo-Systems  model 
1210W-T1.5  probes  (4.5  micron  tungsten)  operated  as  hot  wires  with  a  50% 
overheat.   The  power  and  bridge  for  the  probes  was  the  Thermo-Systems 
model  1054  constant  temperature  anemometer,  the  output  of  which  was 
processed  by  a  Thermo-Systems  model  1057  signal  conditioner  to  remove 
the  dc  level.   The  dc  level  could  cause  the  tape  recorder  to  saturate 
when  large  fluctuation  signals  occured.   The  signal  conditioner  also 
contained  band  pass  filters,  which  were  set  for  a  range  of  0.2  Hz  to 
1  kHz,  reducing  electronic  noise.   The  50%  overheat  used  was  high  enough 
so  that  the  system  response  to  temperature  fluctuations  was  well  below 
the  electronic  noise  level.   The  hot  wire  sensors  were  oriented  with  the 
axis  of  the  wire  in  the  vertical  direction.   With  this  orientation  the 
sensor  was  insensitive  to  the  vertical  component  of  the  wind,  so  that 
only  horizontal  fluctuations  were  recorded. 

The  sensors  for  air  temperature  fluctuations  were  Thermo-Systems 
model  1210W-P.8  (2.5  micron  platinum)  operated  as  cold  wires.   An  ac 
resistance  bridge,  Sylvania  model  140  thermosonde,  which  operated  at 
3  kHz,  was  adapted  for  shipboard  use  and  acted  as  the  detector.   With 
this  unit  data  could  be  obtained  up  to  1  kHz.   It  is  necessary  to  use 
an  ac  bridge  as  the  detector  so  that  high  sensitivity  could  be  obtained 
with  a  very  low  current  through  the  sensor.   If  the  temperature  of  the 
sensor  was  elevated  due  to  bridge  current  the  system  would  sense  wind 
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speed  fluctuations  as  well  as  temperature  fluctuations.   The  bridge  had 
a  long  time  constant,  automatic  balance  circuit  so  that  slow  drifts  of 
ambient  temperature  were  compensated. 

The  thermosonde  configuration  was  such  that  external  resistances 
would  be  used  in  both  arms  of  the  bridge.   This  allowed  one  to  use  a 
single  sensor  balanced  by  a  precision  decade  resistance,  or  to  use  two 
cold  wire  sensors  in  the  bridge.   Thus,  two  methods  were  available  for 
evaluating  temperature  fluctuations:   (1)  A  single  sensor,  the  signal 
being  analyzed  by  a  spectrum  analyser-T'  mode,  and  (2)  two  sensors  spaced 
at  a  fixed  separation-AT  mode,the  resulting  signal  being  processed  by 
an  RMS  module,  with  the  average  either  recorded  on  a  strip  chart  or  by 
the  MIDAS  system. 

Figure  5  shows  that  the  wind  speed  fluctuation  signals  were  also 
processed  by  an  RMS  module  after  the  signal  had  passed  through  a  band 
pass  filter.   The  band  pass  filter  operated  in  the  time  domain  and  per- 
formed a  function  equivalent  to  the  spatial  filtering  accomplished  by 
the  probe  separation  utilized  in  the  AT  mode  for  temperature  fluctuations. 
Spatial  filtering  could  not  be  used  for  wind  speed  fluctuations  because 
it  was  not  possible  to  operate  two  hot  wires  by  a  single  bridge.   Using 
two  bridges  and  a  different  circuit  introduced  too  much  error  into  the 
process . 

A  final  arrangement  for  evaluating  fluctuation  signals  led  to  the 
differentiation  of  the  bridge  output  prior  to  RMS  analysis.   Dif- 
ferentiation necessarily  emphasizes  high  frequency  components  so  that 
high  frequency  noise  becomes  a  problem.   To  alleviate  this  problem  a 
1  kHz  active  filter  with  a  48  dB  rolloff  was  used  but  noise  problems 
were  still  severe  enough  so  that  no  satisfactory  results  were  obtained 
with  this  method. 
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Methods  for  analyzing  data  gathered  utilizing  the  various  experi- 
mental techniques  will  be  described  in  Section  VI. 

Humidity  fluctuations  were  detected  with  a  Lyman-a  cell.   Difficul- 
ties were  experienced  in  keeping  this  system  in  operation  in  the  rigors 
of  a  marine  environment  and  no  results  will  be  reported  here. 

C.   Auxiliary  Equipment 

An   Aerovironment  model  300  Acoustic  Radar  was  installed  on 
the  R/V  Acania  in  order  to  monitor  the  height  of  the  temperature  inver- 
sion.  The  sounder  transmitted  a  1600  Hz  pulse  of  sound  which  was 
scattered  by  small  scale  ambient  temperature  inhomogeneities .   Using  a 
facsimile  recording  device,  the  scattered  return  soundwaves  were  then 
detected  to  produce  the  observed  signal.   The  large  temperature  gradient 
and  strong  turbulence  at  the  base  of  an  inversion  produced  a  strong 
reflected  signal,  so  that  the  height  of  the  inversion  could  be  readily 
determined.   It  was  possible  to  measure  the  strength  of  thermal  turbu- 
lence with  the  acoustic  sounder,  but  noise  and  calibration  problems  made 
such  measurements  impossible  with  our  installation. 

A  Meteorology  Research,  Inc.  model  1580A  Fog  Visiometer  was  installed 
on  the  bridge  of  the  R/V  Acania.   The  instrument  measured  the  amount  of 
light  scattered  from  particulates  and  had  a  usable  visibility  range  of 
8  m  to  13  km.   The  unit  could  be  rotated  so  that  the  relative  wind  was 
incident  directly  on  the  open  side  of  the  scattering  region.   Rotation 
was  necessary  when  the  sun  was  low  on  the  horizon  in  order  to  keep  sun- 
light out  of  the  receiving  optics. 
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IV.   Calibration 

A.   Mean  System 

1.   Temperature 

The  mean  temperature  sensors  (HP  model  2833  Quartz)  were  cal- 
ibrated using  standard  laboratory  thermometric  techniques.   Since  the 
system  was  used  to  measure  the  gradient  of  a  temperature  profile,  the 
relative  consistency  of  the  sensors  was  considered  more  important  than 
the  absolute  temperature  calibration.   All  five  sensors  were  greased  and 
screwed  into  a  3  inch  aluminum  cylinder  fitted  with  a  center-mounted 
thermocouple.   The  aluminum  cylinder  was  0-ring  fitted  into  a  copper 
pipe  long  enough  (24  inches)  to  permit  coiling  of  the  majority  of  the 
armored  sensor  cable  in  the  temperature  controlled  environment.   The 
pipe  and  cylinder  unit  were  immersed  in  a  15  liter  glass  dewar  filled 
with  a  mixture  of  water  and  ice.   The  cylinder  and  pipe  combination 
allowed  the  sensors  and  cables  to  be  surrounded  by  the  bath  without 
actually  getting  wet.   The  large  mass  and  excellent  thermal  conductivity 
of  the  copper  and  aluminum  eliminated  transient  fluctuations  and  reduced 
inter-sensor  temperature  differences  to  a  few  thousandths  of  a  centi- 
grade degree. 

The  normal  calibration  procedure  was  to  cool  the  sensors  to  about 
15  C  while  monitoring  their  indicated  temperatures  with  printout  from 
MIDAS  at  one  minute  intervals.   Once  equilibrium  was  reached,  the  sensor 
oscillators  were  adjusted  so  that  (if  possible)  all  sensors  read  within 
.01  C  of  each  other.   The  temperature  was  then  varied  in  2-3  C  incre- 
ments by  adding  ice  or  warm  water  and  noting  the  MIDAS  printout  values 
when  a  new  equilibrium  was  reached.   Initially,  the  temperatures  were 
compared  to  the  thermocouple  values  (data  set  1/7/76  in  Table  I)  but 
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TABLE  I   Mean  Temperature  Calibrations 
T,  °C 


TC 

0 

1 

2 

3 

4 

.01 

.02 

.00 

.02 

-.01 

.02 

5.25 

4.94 

4.93 

4.94 

4.92 

4.94 

1/7/76        7.50 

7.24 

7.24 

7.25 

7.23 

7.24 

9.45 

9.21 

9.22 

9.22 

9.21 

9.21 

11.95 

11.64 

11.66 

11.65 

11.65 

11.65 

13.35 

13.05 

13.08 

13.06 

13.06 

13.06 

14.85 

14.53 

14.56 

14.54 

14.55 

14.54 

23.25 

22.86 

22.89 

22.87 

22.88 

22.87 

8.61 

8.60 

8.62 

8.61 

8.62 

2/9/77 

13.14 

13.14 

13.15 

13.15 

13.16 

20.54 

20.50 

20.51 

20.52 

20.53 

this  was  later  deemed  unnecessary.   One  problem  with  the  system  developed 
due  to  using  quartz  sensors  with  different  calibration  properties.   In 
data  set  1/7/76  of  Table  I,  sensors  1  and  3  are  from  a  different  batch 
than  0,  2  and  4.   This  led  to  increasing  disagreement  at  the  extremes 
of  the  calibration  range.   This  problem  was  alleviated  by  ordering  more 
sensors  and  specifying  characteristics  to  match  sensors  0,  2  and  4  (data 
set  2/9/77  of  Table  I).   With  this  improvement,  the  accuracy  of  a  single 
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sensor  (relative  to  the  average  of  all  five)  was  no  worse  than  ±  .02  C 

over  a  reasonable  operating  range.   Note  that  this  is  the  accuracy 

of  the  temperature  sensor  only  and  does  not  include  errors  introduced  by 

the  aspirators.   The  relative  aspirator  temperature  error  depends  most 

strongly  on  the  amount  and  uniformity  of  solar  heating;  this  error  is 

quite  small  at  night  or  under  heavy  clouds  but  may  be  as  great  as 

±  .1  C  in  bright  sunlight  with  light  winds. 

2.   Humidity 

To  determine  their  accuracy  and  reproducibility,  detailed  exam- 
ination in  laboratory  conditions  has  been  made  of  the  calibration  char- 
acteristics of  the  Hygrodynamics  Digital  Hygrosensor.   The  calibrations 
were  based  on  ambient  humidity  variations  (from  24%  RH  to  50%  RH  over  a 
one  month  period) ,  as  measured  with  a  Bendix  Psychron  Wet-Dry  Bulb 
Psychrometer,  and  controlled  humidity  (from  40%  RH  to  99%  RH) ,  as  meas- 
ured by  a  wet-dry  thermocouple  in  a  2  x  3  x  3  ft  fog  chamber  adapted  for 
these  experiments.   The  air  circulated  through  the  chamber  could  be  pre- 
humidified  in  a  column  of  wet  glass  beads  giving  nominal  humidity  of 
about  80%.   High  humidities  (up  to  99%  RH)  could  be  obtained  by  intro- 
ducing varying  amounts  of  water  in  the  form  of  a  fine  mist  into  the 
chamber.   In  general,  the  humidity  in  the  chamber  could  be  stabilized 
for  periods  of  several  hours  and  was  uniform  throughout  the  chamber  to 
within  1%  RH. 

The  hygrosensor  was  a  Dunmore-type  lithium  chloride  sensor  whose 
resistance  varies  in  proportion  to  the  relative  humidity  to  which  it 
is  exposed.   It  is  a  slow-response  sensor  and  is  generally  utilized  in 
mean  system  measurements.   The  standard  procedure  of  calibrating  these 
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sensors  in  a  small  chamber  with  a  saturated  salt  solution  of  known  vapor 
pressure  was  not  .satisfactory  for  a  program  of  calibration  of  eight  sen- 
sors.  This  was  primarily  due  to  a  lack  of  confidence  in  our  ability  to 
reproduce  humidity  conditions  in  consecutive  sensor  calibrations.   Attempts 
to  calibrate  all  eight  sensors  simultaneously  in  a  large  volume  chamber 
were  complicated  by  temperature  drifts  that  made  it  difficult  to  assure 
equilibrium  conditions  necessary  to  the  accuracy  of  the  technique.   As 
a  result,  it  was  decided  to  use  the  humidity  controlled  chamber  previously 
described.   The  quantitative  results  of  a  series  of  calibrations  are  given 
in  Table  II;  in  general,  the  eight-sensor  average  agreed  with  the  psychro- 
meter  standard  to  -0.4  ±  2.9%  RH  for  measurements  taken  with  equilibrium 
times  of  order  one  hour.   The  individual  total  humidity  accuracy  of  the 
sensor  was  of  less  importance  than  the  relative  consistency  within  the 
group  of  sensors,  since  they  were  used  in  a  four-level  system  to  determine 
humidity  gradients.   In  this  respect,  they  could  be  calibrated  in  the 
laboratory  to  ±  1%  RH,  subject  to  assurances  of  sufficient  time  for 
response. 

The  suitability  of  the  hygrosensors  for  shipboard  measurement  of 
humidity  gradients  was  reduced  by  two  inherent  properties  of  the  device. 
An  aspirated,  teflon-coated  sensor  has  a  time  constant  of  about  30  minutes, 
but  this  varies  from  sensor  to  sensor.   Therefore,  the  humidity  differ- 
ence read  between  two  sensors  will  be  unreliable  if  the  humidity  is 
changing  faster  than  about  10%  RH  per  hour.   A  more  serious  deficiency 
appears  if  the  sensors  are  exposed  to  humidities  above  95%  RH.   At  those 
high  humidities,  the  sensors  become  sluggish  and  exposure  to  99%  RH  may 
cause  the  sensors  to  be  useless  for  as  long  as  several  days.   Further 
uncertainties  are  introduced  by  the  use  of  long  cables  in  the  shipboard 
system  and  by  the  accumulation  of  sea  salt  on  or  near  the  sensors. 
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TABLE  II   Hygrodynamics  Digital  II  Calibrations 


Wet -Dry 

Ave 

Bulb 

Hygro 

%RH 

%RH 

1 

2 

3 

4 

5 

6 

7 

8 

33. 

31.6 

1.4 

-1.4 

-0.6 

1.5 

-0.2 

-1.6 

-0.9 

1.8 

36. 

37.2 

-0.5 

-0.8 

-0.3 

1.6 

0.5 

-1.4 

-0.1 

1.2 

44.  - 

41.2 

-0.8 

-0.4 

-0.3 

0.6 

0.9 

0.4 

-0.1 

0.3 

NA 

41.5 

-0.5 

-0.4 

-0.3 

0.5 

0.9 

-0.4 

-0.1 

-0.1 

NA 

42.6 

-1.2 

-0.5 

-0.5 

-0.9 

1.8 

-0.4 

-0.2 

0.1 

NA 

65.8 

-0.2 

-1.4 

-2.5 

0.8 

1.8 

-0.3 

1.3 

0.3 

74.8 

73.6 

0.2 

0.2 

-1.3 

0.5 

-0.2 

0.4 

-0.4 

0.8 

88.3 

85.0 

-0.7 

0.9 

-0.4 

0.1 

-0.6 

0.0 

0.4 

0.2 

82.6 

86.7 

-1.4 

-0.6 

-0.5 

0.1 

-0.2 

0.3 

0.7 

0.3 

82.6 

87.0 

-1.0 

0.5 

-0.7 

0.0 

0.3 

0.7 

0.0 

0.2 

85.0 

83.0 

-0.4 

1.1 

-0.1 

0.1 

-1.8 

1.1 

-0.6 

-0.1 

93.9 

89.6 

0.1 

0.8 

-0.2 

0.3 

-2.2 

2.1 

-1.5 

0.4 

94.4 

92.3 

-0.6 

1.1 

-0.5 

-0.3 

-0.4 

1.2 

-0.5 

0.4 

94.4 

93.6 

-0.8 

1.0 

-0.8 

-0.3 

0.4 

-0.2 

0.8 

-0.4 

Ave  Devi at. 

-0.5 

-0.0 

-0.6 

0.3 

0.1 

0.1 

-0.1 

0.4 

a 

0.7 

0.9 

0.6 

0.7 

1.1 

1.0 

0.7 

0.6 

The  average  deviation  of  the  wet -dry  bulb  from  the  Ave  Hygro  reading 
0.4%  RH  with  a  standard  deviation  of  2.9%  RH. 


A-16 


3.   Velocity 

The  mean  wind  velocity  sensors  (Thornthwaite  cup  anemometers) 
were  calibrated  in  a  wind  tunnel  with  a  nominal  one  meter  square  test 
section.   The  cups  were  calibrated  four  at  a  time  with  a  symmetrical 
placement  to  reduce  relative  velocity  differences.   A  summary  of  the 
calibration  is  shown  in  Figure  6  expressed  as  individual  cup  velocity 
deviations  from  the  wind  tunnel  velocity.   From  this  data  it  can  be 
concluded  that,  for  gradient  purposes,  the  cup  anemometers  were  accurate 
to  about  ±  4%.   The  deviation  of  the  calibration  curve  at  low  wind 
speeds  is  due  to  errors  in  the  determination  of  wind  tunnel  velocity 
values. 

B.   Fluctuation  System 

The  essence  of  the  fluctuation  calibrations  is  contained  in 
establishing  the  relation  between  the  instrumental  voltage  variation,  V', 
about  mean  voltage,  V.  with  respect  to  the  atmospheric  fluctuation,  X', 
about  its  mean  value  X,  i.e. 

X'   =   F  (  V,  X  )  V  (1) 

The  calibration  factor  F  (  V,  X  )  may  be  a  function  of  individual 
sensors  as  well  as  the  instrument  being  used.   This  process  is  usually 
accomplished  by  measuring  the  dependence  of  X  on  V  and  relating  it  to 
F  (  V,  X  )  through  the  derivative, 

F  C  v,  x  )  =  f iX®. ]  ~l  (2) 
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1.   Temperature 

The  appropriate  temperature  fluctuation  expression  is; 


T'   =   (  a  G  R  )_1  V  (3) 


where  a  is  the  temperature-resistivity  coefficient  of  the  platinum  sensor 

wire,  G  is  the  gain  factor  of  the  thermosonde  (set  at  G  =  20  volt/ft) ,  and 

R  is  the  resistance  of  the  sensor.   In  a  previous  work  (Schacher  and 

Fairall,  1976)  it  was  found  that  a  =  .0036  °C~ '  with  little  variation 

from  sensor  to  sensor.   The  frequency  response  of  the  thermosonde  units 

was  about  1  kHz. 

2.   Humidity 

The  Lyman-alpha  sensor  measures  humidity  as  a  function  of  the 

o 
absorption  of  ultraviolet  light  at  the  1215A  Lyman-alpha  transition  of 

hydrogen  in  water  vapor.   It  consists  simply  of  a  UV  source  tube  and 
detector,  separated  by  an  absorbing  air  gap. 

The  detector  voltage  is  proportional  to  the  total  light  trans- 
mitted and  is  related  to  the  absolute  humidity,  q  (in  mb) ,  by  the 
equation 

V  =  V  exp  (-Xq)  (4) 

The  Lyman-alpha   is  a  fast  response  device  and   is   specifically  used  for 
measurement   of   humidity   fluctuations,    q' 

q'=       4     ^  (5) 

A      V 

The   quantity  V      in  Eq.    (4)    is   the  dry-air  voltage  which  can   be   found   by 
bathing   the  detector   in  dehumidified   air.      Since  water   soluble  windows 
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(LiF  or  MgF~)  are  required  for  transmission  in  the  ultraviolet,  V   is  a 
highly  variable  quantity,  a  property  which  makes  the  Lyman-alpha  a  poor 
device  for  measurement  of  total  humidity,  q.  However,  Eq.  (5)  shows  that 
the  fluctuating  part,  q',  depends  only  on  X,  which  is  a  property  of  water 
vapor  and  the  gap  setting  only  and  is  independent  of  the  window  condition 
so  that  q'  can  be  reliably  measured.   The  Lyman-alpha  was  calibrated  in 
dry  air,  ambient  air,  and  in  the  humidity  chamber  and  it  was  found  that 
X     ■   .225  mb  "  for  source  =  50  Ua  and  a  1  cm  air  gap  (Fig.  7) . 

Feasibility  of  the  Lyman-alpha  sensor  for  shipboard  humidity  fluc- 
tuation measurements  is  subject  to  considerations  of  the  survivability 
of  the  windows  in  the  ocean  environment  of  high  humidity,  sea  spray,  and 
rain.   In  the  laboratory  humidity  chamber,  V  was  found  to  be  reduced  by 
a  factor  of  four  after  a  three-day  exposure  to  85%  humidity.   Consequently, 
it  was  necessary  to  provide  physical  protection  for  the  sensors  during 
routine  operations. 

3.   Velocity 

The  hot  wires  (or  hot  films)  were  calibrated  using  a  TSI  1125 
calibrator.   It  was  important  to  use  the  actual  cables  from  the  shipboard 
system  (or  their  equivalent)  to  prevent  errors  due  to  voltage  drops  in 
the  cable.   The  appropriate  mean  dependence  function  is 


V   (U)   =   (V  2  +  B  /U  )h  (6) 

o 

2 
where  V  and  B  are  obtained  for  each  sensor  at  a  specific  overheat.   The 
o 

derivative  of  (6)  yields  the  fluctuation  calibration  factor 

0.    =     i-L/J-      v  (7) 
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The  frequency  response  of  the  velocity  bridge  was  about  20  kHz.   The 
sensor  frequency  response  depends  on  the  sensor  diameter  and  overheat, 
but  was  typically  greater  than  1  kHz  (Fairall  and  Schacher,  1977). 
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V.   Data  Acquisition 

Three  criteria  applied  when  deciding  the  methods  for  acquiring  data: 
1)   data  was  to  be  taken  for  a  wide  range  of  local  meteorological  con- 
ditions,  2)   a  given  condition  should  be  tested  several  times,  preferably 
on  different  days  and/or  locations  to  test  the  statistical  validity  of  the 
data,  and,  most  importantly,  3)   shipboard  conditions  must  be  such  that 
valid  data  can  be  obtained. 

The  most  important  consideration  to  insure  validity  of  the  data  was 
the  airflow  over  the  ship.   Since  all  of  the  sensors  were  located  close 
to  the  bow  of  the  ship  the  relative  wind  had  to  be  from  the  bow  or  the 
ship  could  severely   influence  the  data  obtained.   It  was  found  that 
data  could  be  taken  when  the  wind  was  within  30  of  the  bow.   When  the  wind 
was  more  than  45  from  alignment  with  the  bow  the  wind  profiles  were 
noticeably  affected.   The  position  of  the  sun  was  a  second  consideration. 
When  the  sun  was  low  on  the  horizon  and  behind  the  ship  the  aspirators 
for  levels  1  and  2  were  in  shadow  and  those  for  levels  3  and  4  were 
in  sunlight.   Under  this  condition  the  temperatures  at  the  upper  levels 
were  elevated  enough  to  make  the  obtained  temperature  profiles  invalid. 

Since  data  was  acquired  on  a  fluctuating  system,  a  time  average  had 
to  be    performed  to  obtain  final  parameters.   Averaging  over  10  minute 
and  20  minute  time  spans  was  accomplished  depending  on  the  situation. 
Thus,  data  had  to  be  acquired  in  such  a  way  that  the  conditions  remained 
constant  over  a  time  at  least  as  long  as  the  planned  average  time.   This 
means  that  there  had  to  be  close  coordination  between  scientific  personnel 
and  ship  crew  so  that  no  course  changes  occurred  during  a  data  taking  run. 
Of  course,  changes  in  local  conditions  also  affected  the  data,  such  as 
moving  through  the  edge  of  a  fog  bank.   Even  seemingly  inconsequential 
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occurrences  such  as  moving  into  the  lee  of  an  island  when  the  ship  was 
several  miles  offshore  could  affect  the  data.   A  running  log  of  all 
ship  maneuvers  and  position,  and  of  local  conditions  was  maintained  and 
no  averages  were  performed  for  which  it  was  suspected  that  any  condition 
changed  during  the  time  period  under  consideration. 

Several  tests  were  made  and  it  was  found  that  data  could  be  taken 
both  with  the  ship  at  rest  and  underway.   However,  some  of  the  subsequent 
calculations  then  required  that  the  true  wind  over  the  surface  of  the 
water  be  known  accurately.   Uncertainties  as  to  the  ship's  actual  speed 
make  calculation  of  true  wind  from  the  relative  wind  and  ship  heading 
and  speed  too  inaccurate.   Thus  the  ship  was  stopped  frequently  in  order 
to  obtain  the  true  wind. 
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VI.   Data  Analysis 

The  analysis  procedures  changed  considerably  as  the  instrumentation 
evolved  towards  increased  automation.   Averaging  of  mean  data,  initially, 
a  painstaking  hand  process,  is  now  performed  by  MIDAS  at  specified  time 
intervals  (typically  10,  20  or  30  minutes) .   Rather  than  describing  all 
techniques  that  have  been  used,  discussion  is  limited  to  those  of  more 
recent  application. 

A.   Mean  Profile  Analysis 

The  mean  profiles  and  gradients  were  established  under  the  assumption 
of  the  logarithmic  height  dependence. 

X  =  A  log  Z  +  B  (8) 

The  gradient  at  the  height  Z  is 

f  -  f  <9> 

dz     Z 

In  the  case  of  temperature,  one  is  interested  in  the  potential  temperature, 

6, 

9"  =  T  +   .0098  Z  (10) 

and  the  virtual  potential  temperature,  9  , 

"9  =  "9  +   .61  q  T  (11) 

v 

where  q  is  the  specific  humidity  in  grams  of  water  vapor  per  gram  of  dry 
air.   The  actual  fits  of  the  data  were  done  either  by  the  standard  least 
squares  method  or  by  subjective  ("eyeball")  graphical  techniques.   The 
presence  of  individual  erroneous  values  due  to  printer  errors,  micropro- 
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cessor  errors,  and  instrumental  glitches  required  the  data  to  be  subjec- 
tively edited  even  in  the  least  squares  analysis.   In  either  case,  one  is 
required  to  identify  and  reject  "bad"  values  -  a  delicate  and  unsatisfying 
process.   Figure  8  shows  examples  of  profiles  in  temperature,  humidity  and 
velocity. 

B.   Turbulence  Analysis 

There  are  four  different  methods  available  to  calculate  turbulence 
parameters  from  the  fluctuations  in  the  atmospheric  variables:   spectral, 
derivative  (dissipation),  difference,  and  RMS.   The  spectral  method  is 
based  upon  the  assumption  of  "local  isotropy"  and  the  Kolmokorov  -5/3 
slope  of  the  one-dimensional  power  spectral  density,  $    (k) , 

2   -5/3 
4>  (k)   =   .25  C   k  (12) 

X  X 

where  C   is  the  structure  function  for  the  parameter  x  and  k  is  the  wave 
number.   Using  Taylor's  "frozen  turbulence"  hypothesis  (k  =  27Tf/U),  C 
is  found  by  performing  a  fourier  spectrum  analysis  of  a  signal  in  the  fre- 
quency domain  (f ) , 

-2/3 
f  4>x  (f)   *  k  4>x  (k)  =  .25  Cx2  /  ££)  (13) 


Therefore, 


2/3 

(14) 


<*    '-"(f)       [*5/3*x«>] 


which  is  expected  to  apply  in  the  frequency  range  from  about  0.1  to  100  Hz, 
(This  is  the  inertial  subrange  of  turbulence.)   The  structure  function  can 
also  be  found  by  measuring  the  variance  of  the  difference  in  the  variable 
X  at  two  points  separated  by  a  known  distance,  d, 


Cx2  =   [X  (r)  -  X  (r  +  d)]2  d"2/3      (15) 
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2 
Defined  in  this  manner,  C   is  independent  of  d  in  the  inert ial  subrange. 

The  RMS  method  is  based  upon  measuring  the  variance  of  the  signal  fluc- 

2 
tuations  between  selected  frequency  limits  and  relating  this  to  C  "  through 

X 


Eq.     (13) 


k 
(    U 


<f>      (k)    d   k     =     X'2     =      (X'        )2  (16) 


rms 


J    kt 
Again  assuming  k     =      2irf/U,    then 


,  s     ,  _    .2/3         (X'         )2 

r   2  o    /  2tt   I  rms 


D    /        (f^2/3   -   fu-2/3)  (17) 


where  f  and  f»  are  the  upper  and  lower  frequency  limits,  and  X' 


rms 


represents  the  measured  RMS  fluctuations  of  the  filtered  signal.   Since 

the  rate  of  velocity  turbulence  dissipation,  £,  is  traditionally  used, 

2 
rather  than  C   ,  one  should  use  the  relation 

u 

2  2/3 

C    =  2.0  e  '  (18) 

to  calculate  E.   Due  to  the  extreme  sensitivity  to  high  frequency  noise,  it 

was  found  that  the  derivative  method  was  unsatisfactory  for  shipboard  data. 

2    2 
The  spectral  method  was  used  to  calculate  CL  ,  C   and  £  (data  shown 

T    q 

in  Figure  9) .   This  method  has  the  advantage  that  the  quality  of  the  data 
is  immediately  obvious  and  noise  spikes  (60,  120,  130  Hz)  can  be  ignored. 
The  primary  disadvantage  of  the  spectral  method  is  the  time  required  to 
analyze  four  levels  of  data  with  two  or  three  signals  for  each  level.   Al- 
though it  requires  two  microthermal  sensors  per  level,  the  difference  method 

2 
has  proven  to  be  an  ideal  technique  for  measuring  C   because  the  RMS  temp- 
erature difference  fluctuations  (Eq.15)  can  be  automatically  averaged  using 
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the  ADC  capability  of  MIDAS.   Since  this  is  a  temperature  difference 
measurement,  it  is  less  subject  to  external  sources  or  error  (such  as 
humitity  -  salt  induced  temperature  spikes  (Friehe,  1977)).   The  RMS 
technique  has  proven  to  be  quite  satisfactory  for  measuring  £,  since  it 

too  lends  itself  to  automatic  averaging  on  MIDAS.   Comparison  of  spectral 

1/3 
and  RMS  measurements  of  £,  through  the  quantity  U^  =   (K  Z  e)  '  , 

(Figure  10)  are  quite  consistent.   UA  is  the  friction  velocity,  K  is 

Von  Karmon's  constant  (.35)  and  Z  is  the  height  above  the  surface.   This 

expression  is  valid  for  near  neutral  stability  where  most  of  the  data 

was  taken. 
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VII.   Experiments  Performed 

Observational  experiments  were  conducted  aboard  the  R/V  Acania 
(Fig.  1)  during  cruises  to  San  Nicholas  Island  (SNI) ,  in  the  vicinity 
of  SNI,  in  Monterey  Bay,  and  on  an  open  ocean  cruise  in  the  fall  of  1976. 
Dates  and  locations  of  the  experiments  which  yielded  results  presented 
here  are  listed  in  Table  III. 

The  route  taken  to  and  from  SNI  appears  in  Figure  11  for  that  portion 
of  the  cruise  for  which  data  was  taken.   The  ship's  anchorage,  position, 
and  track  in  the  vicinity  of  the  northwest  tip  of  SNI  appear  in  Figure  12. 
For  Monterey  Bay  results,  all  data  were  collected  while  the  R/V  Acania  was 
anchored.   The  general  anchorage  locations  which  were  used  in  Monterey 
Bay  appear  in  Figure  13.   The  anchorage  locations  were  changed  to  accom- 
modate coincident  ship  to  shore  or  shore  to  shore  optical  propagation 
experiments.   Wind  directions  during  the  Monterey  Bay  experiments  were 
generally  from  the  north  to  northwest  so  the  location  of  anchorage  is  not 
expected  to  be  a  factor  in  the  results  and  they  are  not  distinguished  with 
regard  to  location.   Details  pertaining  to  the  1976  open  ocean  cruise  appear 
in  a  separate  report,  Fairall,  et  al.  (1978). 

Changes  were  made  in  the  shipboard  measurement  systems  and  sensor 
mounting  arrangements  during  the  sequence  of  experiments  associated  with 
these  results.   In  general,  mounting  arrangements  were  changed  to  improve 
the  vertical  resolution  of  the  measured  parameters.   During  the  first  SNI 
cruise  (1-4  March  1973) ,  fluctuation  (temperature  and  velocity)  measure- 
ments were  made  at  only  two  levels  and  mean  wind,  temperature  and  humidity 
measurements  were  made  only  at  one  level.   Sea  surface  temperature  was  also 
measured.   During  the  second  SNI  cruise  (19-23  September  1973),  both  fluc- 
tuation and  mean  measurements  were  made  at  four  levels  which  differed 
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slightly  from  the  final  configuration  shown  in  Figure  1.   All  the  Monterey 
Bay  experiments  were  performed  with  the  final  mounting  configurations 
described  above.   Sensors  described  for  fluctuation  and  mean  measurements 
were  used  in  all  experiments  but  the  recording  systems,  particularly  that 
for  the  mean  data,  were  changed. 

Additional  discussions  on  analyses,  instrumentation  and  experimental 
conditions,  along  with  additional  results,  are  contained  in  theses  listed 
in  the  references,  which  utilized  data  from  these  experiments. 
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1973 


1974 


1975 


1976 


1977 


TABLE  III 
SUMMARY  OF  OBSERVATIONAL  EXPERIMENTS  TO  FEBRUARY  1977 


Date 

Location 

Prevailing  Condition 

1-5  Mar 

SNI 

Unstable 

18-21  Sep 

SNI 

Neutral-Stable 

15-19  Oct 

Mtry 

Bay 

Unstable 

26  Nov-5  Dec 

Mtry 

Bay 

Unstable 

15-18  Jan 

Mtry  Bay 

Unstable 

25-26  Feb 

Open 

Ocean 

Stable 

27-29  Feb 

Mtry 

Bay 

Stable 

6-7  Mar 

Mtry 

Bay 

Stable 

13-15  Mar 

Mtry 

Bay 

Stable 

25-28  Mar 

Mtry 

Bay 

Unstable 

17-21  Jun 

SNI 

Neutral 

12-16  Aug 

Mtry 

Bay 

Neutral-Unstable 

16-21  Sep 

Open 

Ocean 

Unstable 

18-22  Nov 

Mtry 

Bay 

Neutral-Unstable 

9-10  Jan 

Mtry  Bay 

Stable 

24-28  Mar 

Mtry 

Bay 

Unstable 

20-24  May 

Mtry  Bay 

Neutral 

25-29  Aug 

Mtry 

Bay 

Neutral 

11-16  Sep 

Open 

Ocean 

Unstable 

9-11  Dec 

Mtry 

Bay 

Unstable 

29  Mar-2  Apr 

Mtry 

Bay 

Unstable 

26-30  Apr 

Mtry 

Bay 

Neutral-Unstable 

26-30  Jul 

Mtry 

Bay 

Neutral 

20  Sep-14  Oct 

Open 

Ocean 

Neutral-Unstable 

15-18  Feb 

Mtry 

Bay 

Neutral 
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VIII.   System  Evaluation 

A.   Mean  System 

Accurate  measurements  of  atmospheric  surface  gradients,  especially 
in  the  marine  boundary  layer,  are  extremely  difficult.   The  severity  of 
the  problem  can  be  illustrated  by  noting  typical  conditions  in  the  marine 
boundary  layer  (Table  IV)  in  terms  of  the  typical  gradients  OX/3Z  at 
Z  =  10  meters) ,  scaling  parameters  (X^)  and  total  difference  from  Z  =  5  to 
Z  ■  15  meters  (AX  =  |x  (15)  -  X  (5)  |).   The  scaling  parameter  is  defined  by 

3X/3Z  =  j|  fx  (19) 

where  f  is  a  stability  corrections  factor  (Businger,  et  al. ,  1971).   To 
measure  X^  to  10%  accuracy  requires  the  single  level  accuracies  indicated 
in  the  table,  compared  to  our  present  accuracies.   The  temperature,  T,  is 
used  in  the  table  rather  than  the  virtual  temperature  because  it  is  the 
variable  actually  measured. 

Although  the  gradients  are  not  required  to  calculate  the  turbulent 
parameters,  they  are  of  utmost  importance  in  evaluating  the  theoretical 
behavior  of  the  boundary  layer  through  specification  of  the  stability 

(Richardson  number)  and  through  normalization  of  the  turbulence  parameter, 

2 
C   (Wyngaard,  et  al.,  1971).   This  lack,  of  gradient  accuracy  led  to  con- 
struction of  a  near  surface  buoy  system  to  exploit  the  Z   dependence  of 
the  gradients.   With  the  present  systems,  the  gradient  measurements  repre- 
sent the  weak  link  in  efforts  to  evaluate  the  turbulent  boundary  layer 
similarity  expressions  in  the  marine  environment. 

B.   Fluctuation  System 

Despite  a  large  array  of  factors  influencing  the  performance  of  the 
turbulence  measurement  system  (noise,  ship  effects,  salt  on  sensors, 
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TABLE  IV 


Typical  Marine  Boundary  Layer  Profile  Parameters 


T,  °C 


U,  m/sec 


q,  gm/kgra 


(3X/32) 


10  met 


-  .020 


.040 


-  .025 


-  .07 


.15 


-  .08 


AX 


.20 


10  % 


5  % 


Present  Accuracy 
Required  Accuracy 


.04 


2  % 


3  % 


1  % 
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calibration  accuracy)  the  values  of  turbulence  parameters  were  quite  con- 

2 
sistent.   Figure  14  is  a  time  history  of  four  levels  of  C   and  Z   showing 

the  expected  height  dependence.   Figure  15  is  a  similar  plot  for  three 
levels  of  z.      The  scatter  inherent  in  the  z   data  is  nicely  illustrated 
by  removing  the  height  dependence  in  a  U^  time  history  (Fig.  15)  of  buoy 
data.  Note  that  the  single  level  scatter  is  much  smaller  than  the  level- 
to-level  disagreements.   The  average  profile  estimate  (based  on  Eq.  19) 
of  UA  is  about  .12  m/sec  whereas  the  average  turbulence  estimate  is  about 
.10  m/sec.   This  is,  in  fact,  good  agreement  and  illustrates  the  superior- 
ity of  the  buoy  system  for  gradient  measurements.   The  compatability  of 
gradient  and  turbulence  measurements  can  be  examined  in  calculations  of 
normalized  heat  flux,  F,  ,  using  the  profile  method 

Fh  =  '  *I  (39/3Z)  (20) 

and  the  turbulence  method 

Fh  -  ±  CT  Z   1/3   F  (Rt)  (21) 

where  K_  and  F(R.)    are  functions  of   stability.      These  methods   have  been 
favorably  compared    (Fairall,    et  al.,    1978)    for   fog   events  during  a  marine 
fog    study  called   CEWCOM-76    (Fig.    16). 

C.      Summary 

It   is  apparent   that   this   line  of   research  is  permeated  with 
instrumental  difficulties.      The  atmospheric   gradients  over   the  ocean  are 
small   enough   so    that   even   state-of-the-art  mean  measurements   give   fairly 
poor  accuracy.      The   turbulence  measurements   are   inherently  well   above 
noise  levels   but   the   sensors   are   subject    to   an  array  of    environmental 
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problems:   salt  and  water  films,  platform  motion,  aerodynamic  flow  dis- 
tortion by  the  ship,  radio  and  radar  interference,  and  even  corrosion. 
Despite  these  problems,  the  system  described  in  this  paper  has  yielded 
valuable  and  unique  data  about  the  open  ocean  boundary  layer.   Those 
theoretical  discrepancies  that  still  exist  will  hopefully  be  resolved 
through  the  continual  evolution  and  improvement  of  instrumentation  and 
measurement  techniques. 
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List  of  Tables 

I .  Mean  Temperature  Calibrations 

II.  Hygrodynamics  Digital  II  Calibrations 

III.  Summary  of  Observational  Experiments  to  February  1977 

* 

IV.  Typical  Marine  Boundary  Layer  Profile  Parameters 
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Figure  Captions 

1.  Naval  Postgraduate  School  Department  of  Oceanography's  R/V  Acania. 

2 .  Buoy  system  for  near  surface  measurements . 

3.  Block  diagram  of  mean  system. 

4.  Schematic  of  aspirated  shelter  and  sensors.   H:  humidity  sensor, 
Q:  quartz  thermometer,  0:  quartz  oscillator,  F:  exhaust  fan, 

R:  radiation  shield. 

5.  Block  diagram  of  fluctuation  system. 

6.  Wind  tunnel  calibration  of  mean  wind  sensors  (cups) ,  sensor  error  vs. 
wind  speed. 

7.  Lyman-alpha  calibration,  output  voltage  vs.  water  vapor  pressure. 

8.  Example  profiles  of  temperature,  humidity  and  wind  speed. 

9.  Example  power  spectra  of  fluctuation  parameters  (<£x(f))  vs.  frequency, 

the  low  frequency  roll-off  is  instrumental  ac  coupling,  the  high  frequency 
roll-off  represents  the  beginning  of  the  dissipation  range:  the  slanted 
lines  represent  Kolmogorov  -5/3  exponential  dependence  a)  temperature , 
b)  humidity,  and  c)  velocity. 

10.  Comparison  of  U  values  from  the  spectral  method  and  the  RMS  method. 

11.  Route  of  R/V  Acania  to  and  from  SNI . 

12.  Positions  and  track  of  R/V  Acania  in  vicinity  of  SNI. 
13  .  Anchorage  locations  in  iMonterey  Bay . 

14.  Turbulence  parameters  vs.  local  time  during  CEWCOM-76  for  Z  =  4.2  m 
(open  circles) ,  Z  =  6.6  m  (solid  circles) ,  Z  =  11.3  m  (triangles) , 
and  Z  =  17.7  m  (x's):  a)  C  2   and  b)  E. 

15.  Friction  velocity  (U#)  vs.  local  time  comparing  three  turbulence  levels 
and  the  gradient  value . 

16.  Calculation  of  sensible  heat  flux,  W  =  C  Fh  where  C  is  the  heat  capacity 
of  air,  using  the  mean  gradients  (open  circles)  and  the  turbulence 
values  (x's) . 
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APPENDIX  B 
MOS  STABILITY  FUNCTIONS 


The  forms  of  the  mean  gradient  functions  (Businger  et  al . ,  1971). 

♦yCe)  -  (i  -  15?)"1/4  c<o 

^(C)   ■   (1   +  4.70  00 

♦T(c)  =  (i  -  i9e)"1/3  5<0 

♦T(5)   =   (1   +  6.40  00 

The  mean  profile  function 

2 
*(5)  =  W(5)  -  31n  *(0  -  21n  (]  Y^)  -  2  tan-1  *(g)  -  tt/2  +  In  (1+^^0 

The  dimensionless  velocity  dissipation  function  (Wyngaard  and  Cote,  1971) 
EjjU)  =  (1  +  .51c2/3)3/2  ?<0 

Ey(0  =  (1  +  2.5C  2/3)3/2  00 

The  dimensionless  temperature  structure  function  parameter  (Wyngaard  et  al , 
1971). 

fTU)  =  4.9(1  -  n  )"2/3  C<0 

fT(5)  =  4.9(1  +  2.4C  2/3)  00 
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APPENDIX  C 
DATA 


C-l 


Table  I-C.  EM  Propagation  data  (bulk  calculations). 


Z*  duct  height    (m) 

Fm  Minimum  Frequency  tripped  (GHz) 

N  Refractivity  (N  Units) 

dN/dZ  Refractivity  gradient  (m~  ) 

7  -7  f\ 

CM  Refractive  index  structure  function  (m"   ) 
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Pm 


dft/dZ 


Cn~2 


jC=.1 
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GMT 


Fm 


du/az 


Cn 


3.5SE-14 

8.10E-14 

-3r55e~14- 
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5..15E-14 
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2  3E-1-1 

38E-14 

3-,-4^fr— H~ 

3.34E-14 

6.93E-14 
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1.37E-13 
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2.26E-13 

6.22S-13 
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1..65E-12 
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L.99E-12 
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3.02E-12 

-2-T6-3-S— 12- 
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3.13E-12 

S^rte-E— 12- 

2.75E-12 

2.59E-12 

53E-12 
1.96E-12 
2.10E-12 
2.52E-12 

2.31E-12 
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1.37E-12 

1.56E-12 
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1 . 51E-12 

i. 34^-12 
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Fp 
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2.35S-12 

2.44E-12 
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2.47E-12 
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2.43E-12 
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■0.247 
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"2T"3"4^=rT2'- 

2. 52E-12 

2.  34 F -12 

-2T?^r-rr 

2.265-12 
2.24E-12 


C-5 


-3  K/d  Z 


2-.  2  55- 12- 

2..34E-12 

2 .  4  QE-12 

1.9  3E-12 

2.00E-12 

1,795-12 

-lrr7  2-E-I2 

1.64E-12 

1.69E-12 
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1.38E-12 

1.405-12 

3rrWr-i-2- 
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Fm 
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;ate 
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Date 
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Frr. 


dft/dz 


Cn    2 
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55 
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•0.17  3 


3 . 425-13 
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58 
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59 
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Table  II-C.  MOS  scaling  parameters  and  EM  propagation  results.  Compari 
son  of  the  bulk  method  (first  line)  with  turbulence  U*  and 
profile  T*,  q*  (second  line). 


Z/L  Monin-Obukhov  stability  parameter,  \ 

U*  Velocity  scaling  parameter  (m/sec) 

T*  Potential  temperature  scaling  parameter  (  C) 

q*  Water  vapor  mixing  ratio  scaling  parameter  (gm/kg) 

Z*  Duct  height  (m) 

Fm  Minimum  frequency  trapped  (GHz) 
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20.3 
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74.2 
13.3 
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-0-0=5 

6:n 

21.9 
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4.6 


7.7 
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7.3 

3.1 


3.9 


12.1 


X_7_ 


4.3 


4.6 


5.4 
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4.2 
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3.3 


79 
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-3.66E-02 

0.50 
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^ 

: 

-2.34E-02 

0.35 

-0.04 

-0.07 

7.0 

19.5 
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# 
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u* 

T* 

q* 

Z* 
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80 

3/0  3 
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-3.58E-02 

0.49 

-0.0  3 

-0.24 

21.5 

3.6 

-2.18E-02 

0.41 

-0.02 

-0.11 

11.2 

9.6 

81 

3/0  3 

800 

-3.36E-02 

0.47 

-0.0  2 

-0.22 

20.2 

4.0 

1.73E-02 

0.40 

0.05 

-0.11 

17.5 

4.9 

82 

3/0  3 

830 

-3.84E-02 

0.43 

-0.01 

-0.27 

23.3 

3.2 

1.61E-02 

0.00 

0.02 

0.01 

0.0 

0.0 

33 

3/0  3 
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-3.31E-02 

0.46 

-0.01 

-0.26 

23.6 

3.1 

3.30E-02 

0.33 

0.05 

0.00 

0.3 

501.3 

34 

3/0  3 
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-4.8GE-Q2 

0.4  1 

-0.02 

-0.26 

21.4 

3.6 

1.08E-01 

0.29 

0.15 

-0.01 

5.3 

25.7 

35 

3/0  3 
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-2.37E-02 

0.  39 

0.0] 

-0.24 

22.6 

3.4 
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0.32 

0.07 

-0.04 

7.5 
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36 

3/0  3 

10  30 

~1.00E~01 

0.36 

-0  .  0  6 

-0.2  5 

17.6 

4.9 

1.48E-01 

0.19 

0.17 

-0.03 

25.0 

2.9 

37 

3/03 

1 10  0 

-4.57E-02 

0  .  3  6 

-0  . 0  0 

-0.25 

21.3 

3.7 

1.41E-01 

0.00 

0.15 

0.01 

2.5 

93.6 

33 

3/03 

1300 

-5.30E-02 

0.46 

-0  .  0  5 

-0.22 

17.9 

4  .3 

-1.04E-01 

0.43 

-0.11 

-0.41 

25.0 

2.7 

39 

3/0  3 

1330 

-3.71E-02 

0.46 

-0.0  2 

-0.22 

19.  5 

4.2 

5.67E-02 

0.38 

0.13 

-0.05 

15.6 

5.9 

90 

3/0  3 

1600 

-1.62E-02 

0.48 

0.01 

-0.23 

23.5 

3.2 

1.16E-02 

0.38 

0.04 

-0.0  9 

13.4 

7.4 

91 

3/0  3 

1630 

-1.48E-02 

0.49 

0.01 

-0.22 

22.8 

3.3 

-6.29E-03 

0.46 

0.02 

-0.17 

19.7 

4.1 

92 

3/0  3 

1700 

-1.85E-02 

0.54 

-0.00 

-0.23 

23.1 

3.  2 

3.32E-03 

0.47 

0.03 

-0.09 

11.6 

9.1 

93 

3/0  3 

1830 

-6.11E-02 

0.52 

-0.10 

-0.18 

14.2 

6.3 

-7.28E-Q2 

0.41 

-0.11 

-0.31 

21.9 

3.5 

94 

3/0  3 

1900 

-5.04E-02 

0.49 

-0.06 

-0.19 

15.7 

5.8 

-4.99E-02 

0.41 

-0.06 

-0.23 

19.3 

4.3 

95 

3/0  3 

1930 

-3.07E-02 

0.53 

-0.04 

-0.18 

16.9 

5.2 

-3.36E-02 

0.48 

-0.04 

-0.24 

21.6 

3.6 

C-15 

# 

96 
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3/0  3 


97 3/0  3 

9  3     3/0  3 

99     3/03 


GMT 
20  00 

20  3  0, 
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2130 


Z/L 

-7.57E-02 
-1.27E-01 

-7. 63E-02 
-5.96E-02 

-7.33E-02 
-6.52E-02 

-8.63E-02 
-6.75E-02 


u* 

0.50 


T* 
-0.11 


0.40 

0.43 
0.37 

0.46 
0.3  3 


-0.18 
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-0.0  5 

-0.0  9 
-0.08 
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0.32 
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-0.44 
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-0.21 

-0.21 

-0-18 
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6.6 
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5.3 


7.1 


4.7 


100 

3/0  3 

2200 

-9.52E-02 

0.38 

-0.03 

.   -  0 .  15 

12.3 

3.3 

-1.0  3E-01 

0.30 

-0.  0  7 

-0.31 

3.3 

101 

3/0  3 

2230 

-3. 96E-02 

0.3  7 

-0  .03 

-0.23 

15.7 

5.3 

-3.35E-02 

0.30 

-0.04 

-0.30 

21.7 

3.6 
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3/0  3 

2  3  00 

-8. 2GE-02 

0.38 

-0,05 

_    -0.2  3 

1 7 .  0 

5.  1 

-7.38E-02 

0.3  4 

-0.04 
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19.1 
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3/03 
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-0.2  0 
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105 
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5.3 
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3/04 
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-1.01E-01 

0.38 
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-0.23 
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3/04 
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-1.03E-01 
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T* 

q* 
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112 

3/0  4 

500 

-1.79E-01 

0.31 

-0.09 

-0.23 

16.9 

5.  2 

-1.  52E-G1 

0.23 

-0.03 

-0.23 

17.3 

5.0 
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3/ J  4 

53  0 

-1. 91E-01 

0.30 

-0.09 

-0.23 

15.6 
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0.23 
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15.4 

6.0 
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3/04 
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-2.10E-01 

0.23 
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-0.28 

16.4 

5.4 
6.2 

-1.16E-01 

0.19 

-0.04 

-0.21 

14.9 

115 

3/04 

630 

-1.62E-01 

0.33 

-0.0  9 

-0.23 

17.3 

5.0 

-1.34E-01 

0.28 

-0.08 

-0.24 

16.1 

5.5 

116 

3/04 
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-1.18E-01 
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-0.29 

19.3 

4.2 

-1.21E-01 

0.27 
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-0.18 

12.7 

8.0 
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3/0  4 
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-1.41E-01 

0.33 
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-0.29 
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-0.11 
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-1.33E-01 
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-0.26 

18.5 

4.5 
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0.47 
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0.35 


0.46 


■0.03 

-0.18 
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0.21 
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23.8 
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20.0 
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0.33 


0.46 
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3.58E-02 
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Fm 
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4.5 

3.9 
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3.1 


3.8 
3.7 
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3.7 
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4.0 
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4.0 
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4.0 


1 

-8.25E-02 

0.33 
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-0.27 
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-7. 35E-02 

0.59 
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3.1 
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-0.0  5 
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17.9 

4.7 
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-7.93E-02 
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0.  36 
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-0.33 


?2    ? 
21.  6 
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3.6 
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3/06 

300 

-2.72E-02 
-9.57E-02 

0.33 
0.56 

-0.02 
-0.18 

-0.27 
-0.33 

24.7 

21.5 

2.9 
3.6 

177 

3/06 
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-3.08E-02 
-1.1SE-01 

0.35 
0.54 
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-0.35 
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21.1 

2.3 

3.7 

. 

-3.28E-02 

0.32 

-0.03 
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Table  IV-C.  Meteorological  data.  The  data  is  given  by  level  number 
(ie.,  U2  is  the  wind  speed  at  level  2). 


U     Relative  wind  speed  (m/sec) 

PHI    Relative  wind  direction  (degrees) 

TS    Sea  surface  temperature  (  C) 


T     Air  temperature  (  C) 
H     Relative  humidity  (%) 
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